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Abstract 

We study the semiclassical time evolution of observables given by matrix val- 
ued pseudo differential operators and construct a decomposition of the Hilbert space 
L^(M'^) (S> C" into a finite number of almost invariant subspaces. For a certain class 
of observables, that is preserved by the time evolution, we prove an Egorov theorem. 
We then associate with each almost invariant subspace of L^(M'^) (^C" a classical sys- 
tem on a product phase space T*IR'^ x O, where O is a compact symplectic manifold 
on which the classical counterpart of the matrix degrees of freedom is represented. 
For the projections of eigenvectors of the quantum Hamiltonian to the almost invari- 
ant subspaces we finally prove quantum ergodicity to hold, if the associated classical 
systems are ergodic. 
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Introduction 



The relation between dynamical properties of a quantum system and its classical limit is a 
central subject in the field of quantum chaos. In this context quantum ergodicity is a well- 
estabhshed concept [|Zel87| , pdV85| , |HMR87| , |Zel96|| . It states for quantisations of ergodic 



classical systems that the phase space lifts of almost all eigenf unctions of the quantum 
Hamiltonian converge in the semiclassical limit to an equidistribution on the level surfaces 
of the classical Hamiltonian. The principal goal of this paper is to establish quantum 
ergodicity in systems whose degrees of freedom can be divided into two classes such that 
they are represented in the Hilbert space L^(]R'^) (g) C". The semiclassical limit shall be 
performed in terms of a parameter h —>■ which is primarily linked to the (translational) 
degrees of freedom that are described by the infinite-dimensional factor L^(R'^). The finite 
dimension n of the other factor is fixed. Examples for systems where this description can 
be applied are relativistic particles with spin 1/2 in slowly varying external fields gov- 
erned by a Dirac-Hamiltonian, or adiabatic situations modelled with a Born-Oppenheimer 
Hamiltonian. 

This setting leads to a representation of quantum mechanical observables as matrix 
valued pseudodifferential operators acting on L^(R°') ® C", whose symbols are suitable 
matrix valued functions on the phase space T*]R'^ = M'' x R'' with an expansion in h. 
In particular, the principal symbol Hq of the selfadjoint quantum Hamiltonian ?i is a 
hermitian matrix valued function on T*R'^. Its spectral resolution requires to introduce 
several classical dynamics on T*M'^, each of them generated by one eigenvalue of Hq. Lifted 
to the quantum level, this structure results in a decomposition of the Hilbert space L^(M'^)(8> 

into almost invariant subspaces with respect to the dynamics generated by the quantum 
Hamiltonian Ti that is directly associated with the spectral resolution of Hq. Recently the 
case of matrix valued operators for certain quantum Hamiltonians with scalar principal 
symbol, such that on the classical side one still has to deal with a single system, has been 
considered in [PGOO| , |BGK01|| . Here we extend this approach to the general setting of 



matrix valued operators where one has to define suitable classical systems corresponding 
to each almost invariant subspace of the Hilbert space. 

So far it appears that the semiclassical limit has only been performed with respect to one 
type of the degrees of freedom. For a complete (semi-) classical description of the quantum 
systems under consideration one would also require the second type of degrees of freedom, 
that are represented by the factor C" of the Hilbert space L^(R'^) ® C", to be transferred 
to a classical level. It turns out, however, that for this purpose no further semiclassical 
parameter is needed and the dimension n of the second factor can be held fixed. Indeed, 
a suitable Stratonovich-Weyl calculus ||Str57|| allows to map the principal symbols (with 
respect to the parameter h) of observables and their dynamics in a one-to-one manner to 
genuinely classical systems associated with the decomposition of L'^{M.'^) ® into almost 
invariant subspaces. On this classical level the hierarchy of the two types of degrees of 
freedom is reflected in the structure of the classical dynamics: these are skew-product 
flows built over the Hamiltonian dynamics generated by the eigenvalues of Hq. 

Apart from classical ergodicity the proof of quantum ergodicity typically requires two 
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essential inputs. The first one is a suitable version of an Egorov theorem ||Ego69|| that allows 
to express the time evolution of quantum observables in the semiclassical limit in terms of 
a classical dynamics of principal symbols. We achieve this in two steps: beginning with 
matrix valued principal symbols, we proceed to a completely classical level by exploiting 
the Stratonovich-Weyl calculus in the form developed in ||FGV90|| . It is in the last step 
where the skew-product flows become relevant. The second input is a Szego-type limit 
formula that relates averaged expectation values of observables to classical phase space 
averages. This can be obtained by a straight-forward generalisation of previous results 
HMR87I, iBGOOj. 



Our main results are the Egorov theorem in section ^ and the quantum ergodicity the- 
orem in section In order to formulate the Egorov theorem we first identify a subalgebra 
in the class of bounded semiclassical pseudodifferential operators that is invariant under 
the time evolution. The operators in this subalgebra have to be block-diagonal with re- 
spect to the projections onto the almost invariant subspaces of L'^{M.'^) ® C". Theorem |3.2| 
then asserts that the (matrix valued) principal symbol of each block is evolved with the 
Hamiltonian flow associated with that block. In addition, it is conjugated with unitary 
transport matrices that describe the time evolution of the matrix degrees of freedom along 
the trajectories of the Hamiltonian flow. 

We next identify the dynamics provided by the transport matrices with a coadjoint 
action of a certain Lie group. Kirillov's method of orbits [|Kir76|] then enables us to connect 
the apparently quantum mechanical dynamics with a genuinely classical dynamics on a 
certain coadjoint orbit O, which is a symplectic manifold. This relation can be constructed 
explicitly with the help of the Stratonovich-Weyl calculus developed in |FUV5U|. As a result 



we obtain that after a Stratonovich-Weyl transform the principal symbol of each block 
of an observable is evolved with a skew-product dynamics on the combined symplectic 
phase space T*R'^ x O. This observation restores the general picture behind Egorov- 
type theorems: in leading semiclassical order the quantum mechanical time evolution is 
determined by classical dynamics. 

The decomposition of L^(M°') (8)C" into almost invariant subspaces and the correspond- 
ing set of classical flows force quantum ergodicity to be concerned with projections of the 
eigenvectors of 7i to the subspaces, since only these are associated with unique classical 
systems. The projected eigenvectors, however, are no longer genuine eigenvectors of H, but 
only provide approximate solutions to the eigenvalue problem and thus yield, after nor- 
malisation, quasimodes (see [ [Laz93|| ). For the latter we prove quantum ergodicity to hold 
in the usual sense. In this context the relevant version of the Egorov theorem introduces 
on the classical side the skew-product flow associated with the given subspace as described 
above. We show that if this flow is ergodic, the phase space lifts of almost all normalised 
projected eigenvectors converge to equidistribution on the product phase space. 



3 



1 Background on matrix valued pseudodifFerential 
operators 

In this section we recall some basic results of pseudodifferential calculus which are well 
known in the context of operators with scalar symbols. They carry over to the case of 
matrix valued symbols by only slight modifications of the results known for operators with 



scalar symbols which can, e.g., be found in [[Rob87|, pS99|l ; for the matrix valued case see 
also 



The quantities we are primarily concerned with are linear and continuous operators 
B : ^(R'') O C" ^ y{R'^) O C" with Schwartz kernels taking values in the n x n 
matrices M„(C). Instead of using a kernel Kj^ G S^'{M.'^ x M'^) M„(C) an operator B can 
alternatively be represented by its (Weyl) symbol B e ^'{T*R'^) ® M„(C) that is related 
to the Schwartz kernel through 

Here h G (0, ^o]; with Hq > 0, serves as a semiclassical parameter and T*R'^ := R'^ x 
M*^ denotes the cotangent bundle of the configuration space R'', i.e., T*R'^ is the phase 
space of the translational degrees of freedom. Below (see section ^ T*R'^ will provide one 
component of a certain combined phase space, which also represents the degrees of freedom 
described by the matrix character of the symbol in terms of points on a suitable symplectic 
manifold. 

According to the Schwartz kernel theorem every continuous linear map B : ^(R*^) (g) 
C" ^'(R"') (g) can be viewed as an operator with kernel of the above form. However, 
operators with kernels in ^'(T^R'^) O M„(C) are too general for many purposes; e.g., they 
can in general not be composed with each other. One therefore has to restrict to smaller 
sets of kernels and hence to smaller classes of symbols. To achieve this we make use of 
order functions which have to fulfill a certain growth property in the 

sense that there are positive constants C, N such that 



m 



(x,0<C{l + {x- yf + (e - V?) m{y, r^) 



for all (x, ^), {y,ri) G T*R . A typical example for such an order function is given by 



mix 



,o = {i+x'+ef, M>o 



This notion allows us to define the symbol classes which we will employ in the subsequent 
discussions (see |pS99|| ). 



Definition 1.1. Let m : T*R'^ (l;Oo) be an order function. Then define the symbol 
class S(m) C C°°(T*R'^) » M„(C) to be the set of S G C~(T*R'^) ® M„(C) such that for 
every {x, C,) G T*R'^ and all a, /3 G Nq there exist constants Ca,i3 > with 

\\d?d^B{x,0\\nxn < Ca,M^,0- (1-2) 
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Here || ■ ||nxn denotes an arbitrary (matrix) norm on M„(C). If in addition the symbol 
B{x,C,',h) depends on the parameter h G (0, /iq]? we say that B G S(m) if B{-,-;h) is 
uniformly bounded in S(m) when h varies in {0,ho\- In particular, for g G M let S'^(m) 
consist of B : T*M'^ x (0, Hq] M„(C) belonging to n-«S(m), i.e., 

\\dp^^B{x,^;h)\Un < C«,^r^m(x,0 

for all a,(3e N^, {x, e T^M'^, and ^ G (0, ho]. 

An asymptotic expansion of 5 G S'^°(m) is defined by a sequence {Sj G S*(m)}jgNQ of 
symbols, where qj decreases monotonically to — oo and 

N 

B-^BjE 8"^+^ (m) 

j=0 

for all N eNo- In this case we write 

oo 
j=0 

In the following we will often use the class S^i(m) of classical symbols, whose elements have 
asymptotic expansions in integer powers of h, 

oo 
j=0 

where Bj G S(m) is independent of h. We also use the notation 

S°°(m) := y S'?(m) and S-°°(m) := p| S'?(m). 

<jeIR geM 

An operator with a kernel of the form ( p, . 1|) and symbol B G S(m) clearly maps both 
^(M*^) and t5^'(]R'^) ® C"' into themselves, whereby according to (IJ.) it acts on 
C"- valued functions ^ G ^(M°') (g) as 

= (op^[i?]^)(x) = I ^^^^^ei(--^)«5(^,e)^(y) d2/de. 

Operators ;B = op '^[5] of this type are called Weyl operators, and symb'^fS] = B denotes 
the Weyl symbol of B. If the Weyl symbol of an operator is a classical symbol with asymp- 
totic expansion B ~ X^jgn,, h~'^~^^Bj, we also call op'^[i?] a semiclassical pseudodifferential 
operator. The leading order term symb^[S] = Bq is then referred to as the principal 
symbol of B, and the subsequent term Bi as the subprincipal symbol. 

The set of Weyl operators with symbols from the classes S(m) is stable under operator 
multiplication, in the sense that the operator product is again a Weyl operator with symbol 
in a certain class: 
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Lemma 1.2. Let mi,m2 be order functions. Then for Bj G S{mj), j = 1,2, the product of 
the corresponding operators Bj = op^[Bj] is again a Weyl operator that can be expressed 
in terms of the symbols -81,-82 



8^62 = op^ [B,] op^ [B2] = op^ [i?l#52], 

where the symbol product (Bi, B2) ^ Bijj^B2 is continuous from S(mi) x S(m2) to S(mim2) 
in the topology generated by the seminorms associated with the estimate (\1.3^. In explicit 
terms the symbol product reads 



y=x ' 



where a{vx, v^; Wx, w^) = Vx ■ — ■ Wx denotes the symplectic two-form on T*M'^. Fur- 
thermore, Bj G S°i(mj) are mapped to Bi^B2 G S°i(mim2) with (classical) asymptotic 
expansion 

(5i#fi2)(x,0~ E ^-]^U^(9x,d^-dy,d,)\ B,,,{x,0B2,Uy,v) 



y=x 



The following result, which in its original version is due to Beals [Pea77|| , is useful in 
situations where one wishes to identify a given operator as a pseudo differential operator. 

Lemma 1.3. Let B{h) : S^{W^) ®C" — S^'{W^) ® C" be a linear and continuous operator 
depending on the semiclassical parameter h G {0,ho]. The following statements are then 
equivalent: 

(i) B{h) = op'^[-B] is a Weyl operator with symbol B G S*^(l). 

(a) For every sequence li{x,^), . . . ,In{x,C,), G N, 0/ linear forms on T*]R'^ the op- 
erator given by the multiple commutator [op^[/Ar], [op^[/Ar_i], ■ • • , [op^[Zi], i3] ■ ■ ■] is 
bounded as an operator on -L^(]R'^) (g) C" and its norm is of the order . 

The direction (i)^(ii) is a simple consequence of the symbolic calculus outhned above. 
For the reverse direction see [|HS88| , |DS99| . 



In the discussions below we will basically encounter two types of (Weyl) operators: 
quantum Hamiltonians H = op^[H] with symbols H G S°i(m) generating the quantum 
mechanical time evolution, and observables B = op^[B]. In typical Hamiltonian 
H is given and one is interested in a suitable algebra of observables that allows to study 
dynamical properties of the quantum system. For this purpose it is often convenient to 
consider bounded operators. In the scalar case it is sufficient to know the boundedness 
of the symbols in order to obtain a bounded Weyl operator. This result, originally going 
back to Calderon and Vaillancourt ||CV71|| , generalises to the context of pseudodifferential 
operators with matrix valued symbols without changes. 
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Proposition 1.4. Let B{h) G S(l), then the Weyl quantisation op^[B{h)] of this symbol 
is continuous on L'^{M.'^) ®C". Furthermore, for h G (0, ^o] there exists an upper bound for 
the operator norm of op^[B{h)]. 

For a proof of this result in tlie context of semiclassical pseudodifferential operators 
(depending on a parameter H) see |[Rob87| , [HS88| , PS99|| . 



A quantum Hamiltonian is required to be (essentially) selfadjoint. Thus, in the case of 
a Weyl operator Ti. = op^[H] one requires the symbol H to take values in the hermitian 
n X n matrices. In order to trace back spectral properties of Ti. to properties of the 
principal symbol Hq we will have to construct (asymptotic) inverses oiTi — z and relate 
them to [Hq — z)~^. In this context an operator B = op^[B] is called elliptic, if its symbol 
B G S(m) is invertible, i.e., if the matrix inverse B~^ exists in S(m~^). In such a case one 
can construct a parametrix Q G S(m~^) which is an asymptotic inverse of B in the sense 
of symbol products. 

To see that such an inverse exists for elliptic operators, consider 

op^[B] op^[5-i] = 1 - nop^[i?], 

with R G S(m). For sufficiently small h the operator 1 — /j,op^[i?] possesses a bounded 
inverse and one can define a (left and right) inverse op'^[-B^^](l — hop^[R])~^ for op^[-B]. 
Furthermore, the Beals characterisation of pseudodifferential operators (Lemma |1.3|) im- 
plies that this inverse is again a bounded pseudodifferential operator, see also ||DS99 |. 



To obtain an asymptotic expansion for the parametrix Q one next defines the opera- 
tor Qn := op^[B-^]{l + hn + --- + h^n^), with 7^ = op^[R], which is equivalent to 
Q = op^lQ] modulo terms of order h^^^ . One can hence write 

g ~ + h{B-^ifR) + h\B-^i^Ri^R) + ■ ■ ■ , (1.3) 

and finally observes: 

Lemma 1.5. Let B G S(m) be elliptic in the sense that B^^{x,C,) exists for all {x,C,) G 
T*]R'^ and is in the class S{m~^). Then there exists a parametrix Q G S{m^^) with an 
asymptotic expansion of the form (]77^ such that 

Bi^Q ~ Q#5 ~ 1. 



From ( |1.3|) one moreover observes that an elliptic operator with classical symbol has a 
parametrix that is a classical symbol. 

Frequently it is very convenient to have a functional calculus of pseudodifferential oper- 
ators available. In some places, e.g., we would like to apply the Helffer-Sjostrand formula, 
which shows that a smooth and compactly supported function / G C^(M) of an essentially 
selfadjoint operator B with symbol in i? G S(m) yields a pseudodifferential operator 

f{B) = -- [ d,f{z){B-z)-' dz, 
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C^(C) is an almost-analytic 



Imzr for all N E Nn 



whose symbol is in S(m~ ) for every N E N. Here / E 
extension of /, i.e., f{z) = f{z) for ^ G M and \df{z)\ < Cn 
In the scalar case these results were shown in |[HS89|| (see also |PS99|| ) and have been 
extended to the matrix valued situation in ||Dim93| , Pim98|| . A criterion that guarantees 
the essential selfadjointness of B is that first its symbol B E S(m) is hermitian and, second, 
that B + i E S{m) is elliptic in the sense described above. If B E S°i(m) one can even write 
down a classical asymptotic expansion for the symbol of the operator f{B) whose principal 
symbol reads /{Bq), where Bq denotes the principal symbol of B, see |[Rob87| , PS99|| . 



2 Semiclassical projections 

We motivate the following construction of semiclassical projection operators by considering 
the time evolution generated by a quantum Hamiltonian Ti, i.e., the Cauchy problem 

ih^^i;{t)=nm (2.1) 

for an essentially selfadjoint operator 7i defined on a dense domain D{l-L) in the Hilbert 
space LF'{W^) ® C". If one introduces the strongly continuous one-parameter group of 



unitary operators lA{t) := exp(— -^Tit), t G M, a solution of (2J.) can be obtained by 
defining ip{t) := W(t)V'o for '^o e D{n). Therefore the time evolution B{t) := U{t)*BU{t) 
of a bounded operator B E B{L'^{W^) ® C") in the Heisenberg picture has to fulfill the 
(Heisenberg) equation of motion 

^B{t)='-[H,B{t)]. (2.2) 

If one assumes B and Ti to be semiclassical pseudodifferential operators with symbols in 
the classes S^j(l) and S° (m), respectively, equation ( ^.2] ) yields in leading semiclassical 
order an equation for the principal symbols: 

^So(t) = -^[i/o, Bo{t)] + 0(n°), n ^ 0. 

If one now requires the time evolution to respect the filtration of the algebra S^(l) := 
Uggz^ci(-'-) then, in particular, the principal symbol -Bo(^) should stay in its class that 
derives from the associated grading S^j(l)/S^]~^(l), q E Z. One thus has to restrict to 
operators whose principal parts Bq commute with the principal symbol Hq of the operator 
7{. This condition is equivalent to a block-diagonal form of Bq, 

I 

5o(a;,0 =I]^M,o(^,05o(a;,0^/.,o(a;,0, (2-3) 

with respect to the projection matrices P^^o • T*]R'^ M„(C), = 1,... onto the 
eigenspaces corresponding to the eigenvalue functions : T*R'^ — R of the hermitian 
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principal symbol matrix Hq : T*R'^ M„(C). Since ( p.3| ) is the semiclassical limit of 
the symbol of the operator op'^[P^,o]'Bop'^[P^,o]) when i3 is a semiclassical Weyl 

operator with symbol B G S^i(l), one can ask how the symbols -P^,o, which are projec- 
tors onto the eigenspaces of Hq in C", are related to projection operators onto (almost) 
invariant subspaces of L^(M'^) (S)C" with respect toTi = op'^[if]. We are hence looking for 
quantisations of symbols G S°i(l), with principal symbols P^^, which are (almost) 
orthogonal projections, i.e., 

V,V, = V^ = V; mod 0{h^) (2.4) 

in the operator norm. Moreover, in order that these operators map to almost invariant 
subspaces of L^(M'^) ® C"" with respect to the time evolution U{t) = exp(— -^Tit) generated 
by Ti, we require them to fulfill 

||[^,^m]||l2{r<^)®C" = modOih'^). (2.5) 



As it will turn out, it is even possible to modify the operators in such a way that they 
satisfy the relation ( |2.4| ) exactly, i.e., not only mod 0{h°°). 

The above requirements lead us to consider (formal) asymptotic expansions for the 
symbols P^, 

oo 
j=0 



which satisfy ( p.4| ) and ( |2.5| ) on a (formal) symbol level: 

p^#p^ ~ p^ ~ p;, (2.6) 

and 

[P„ HU ■■= P.i^H - Hi^P, ~ 0. (2.7) 

The solutions of the above equations will be called semiclassical projections and can be 
constructed by two different methods. The first one is based on solving the recursive 
problem arising from ( |2.6|) and ( p.7|) by employing asymptotic expansions of P^ and H 
in S° (1) and S° (m), respectively, using the symbolic calculus outlined in section |^ and 
finally equating equal powers of the semiclassical parameter h. For this procedure cf. 
tEW96 , BN99 ]. The second method employs the Riesz projection formula in the context of 



pseudodifferential calculus [[HS88| , [NS01|| . In the following we will pursue the latter method. 

To this end we consider the matrix valued hermitian principal symbol Ho G S(m) of 
the operator Ti, and in the following we assume: 

(HO) The (real) eigenvalues A^, /i = 1, . . . ,1, of Hq have constant multiplicities ki, . . . ,ki 
and fulfill the hyperbolicity condition 

\\u{x,^) — \^{x,^)\ > Cm{x,^) for 7^ /i and |x| + |^| > c. 
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This requirement is analogous to a condition imposed in ||Cor82|| on the eigenvalues of the 



symbol of an operator in a strictly hyperbolic system, i.e., where the eigenvalues are non- 
degenerate. In particular, the problem of mode conversion that arises from points where 
multiplicities of eigenvalues change is avoided. Since the eigenvalues are solutions of the 
algebraic equation 

n 

det(i/o(x,0-A) = J]r/,(x,OA'^ = 0, (2.8) 

they are smooth functions on T*M'^. Moreover, since Hq is supposed to be hermitian, 
the eigenvalues are bounded by the matrix norm of Hq. Using the smoothness of the 
eigenvalues and the hyperbolicity condition (HO), one obtains: 

Proposition 2.1. Let H e S[!j(m) he hermitian and let the hyperbolicity condition (HO) 
he fulfilled. Then there exist symhols G S°j(l) with asymptotic expansions 

oo 

P^^Y.^^P,j, fi = l,...,l, (2.9) 



that fulfill the conditions i\2.di ) and \2. In particular, the coefficients P^j, j G No, are 



unique, i.e., the symhols are uniquely determined modulo S~°°(l). 

Furthermore, the corresponding almost projection operators = op^[P^] 
semiclassical resolution of the identity on L^(]R^) ® C'^, 

Pi H \-Vi= idi2(Kd)c>3C" mod 0(h°°). 

Proof. We use the technique of ||HS88| , [NSOlf and consider the Riesz projections 



P^{x,0-=7^[ Q{x,C,z)dz, (2.10) 

where r^(x,^) is a simply closed and positively oriented regular curve in the complex 
plane enclosing the, and only the, eigenvalue X^{x,^) G M of Ho{x,C,). Moreover, Q{x, ^, z) 
denotes a parametrix for H — z, i.e., {H — z)^Q ~ Q^{H — z) ~ 1 that will be constructed 
below. For technical considerations one may choose the contour as r^(x,^) = {X^{x,^) + 
Pf^{x,^) < < 2tt} with < c < < |mini.^^{|A^ — Xu\}- Since Hq is hermitian 
with eigenvalues X,y, v = 1, . . . ,/, one can estimate the matrix norm of {Hq — z)^^ for 
z G T^{x,C) as 

||(i7o(^,0-^)"inxn< ^ 



The condition (HO) then allows to choose p^(a;, ^) > cm(x, ^), so that i^o — 2; is elliptic 
for z G r^. If then h is sufficiently small, also H — z = Hq — z + 0{h) is elliptic and the 
relation 



{H - z)i^{Ho- z) ^ = l-hR 
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enables one to construct a parametrix Q{x,^,z) G S°i(m ^) ioi H — z with asymptotic 
expansion 

oo 

Q{x,^,z) ^^h^Qj{x,^,z) 

j=0 

in the same manner as in ( |1.3|) , see also ||Rob87| , [DS99|| . Plugging this expansion into ( |2.10|) 
one obtains 



^TTl Jr„(x,£) 



(2.111 



/ Q,ix,^,z)dz=:y2WP^,,ix,0 
Jrjx,£) 



by using the Borel construction to sum asymptotic series of symbols. According to the 
properties of the Riesz integral the symbols therefore fulfill ( p.6|) and ( p.7|) . Since these 
equations have unique solutions modulo 0{h°°) |pi)W96|| , the coefficients P^j are unique. 

We now consider more general z ^ C, and by inspecting the above construction notice 
that the parametrix Q{z) is well-defined as long as z has a sufficiently large distance from 
the eigenvalues of Hq. According to equation (|1.3| ) its asymptotic expansion then reads 



Q{z) ~ {Ho - z)-^ + h{Ho - 2)-i#i?(2)#(idcn + hR{z) + h^R{z)i^R{z) + 
Since 

R{z)=^-{l-{H-zmHo-z)-'), 



it follows according to the composition formula of Lemma |1.2| that R{z) contains a factor 
{Hq — z)~^, and therefore the only singularities of Q{z) are caused by the eigenvalues of 
-f^o- Thus, according to the Cauchy formula the expression 



Pi + --- + Pi = ^ ! Q{z) dz 



can be replaced by 



1 

27ri 



Q{z) dz, 
r(r) 



where r(r) is a contour with minimal distance r from the origin in C that encloses all 
eigenvalues of while keeping a sufficient distance from them. The value of the above 
integral does not depend on the particular choice of r(r) so that one can take the limit 
r — i> oo and hence obtains 

lim [ Q{z) dz = lim ^ / (Hq - z)-^ dz = idc- mod 0{h°°). 

□ 



r^oo 2m Jr{r) r^oo 2711 Jr{r) 
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The so constructed symbols yield semiclassical almost projection operators 



op^[P,] 



which according to Proposition O are bounded and obviously satisfy the relations (^) 
and ( |2.5| ). Following |[Nen99|| one can even construct pseudodifferential operators that 
are semiclassically equivalent to in the sense that \\P^ ~ 0{tf°), and which fulfill 



41) exactly. To see this, consider the operator 



\z~l\-- 



which is well-defined since the spectrum of is concentrated near and 1. Thus is an 



orthogonal projector acting on L 



C", with \\[v^,,n]\\ < c\\[Vf,,n]\\ = 0{h^). since 



Vfj, is close to Vfj, in operator norm, Beals' characterisation of pseudodifferential operators 



(see Lemma |1.3|) yields that is again a pseudodifferential operator with symbol in 
the class S''(l). This has already been noticed in ||JNS01|| and follows from the fact that 
{Vn — z)"-^ for 1^; — l| = l/2isa pseudodifferential operator according to the parametrix 



construction of Lemma |1.5| . Having projectors available, one can also construct (pseudo- 
differential) unitary transformations of L^(R'^) ® C" which convert Ti by conjugation in 
an almost block-diagonal form, see ||(Jor83b| , [Ll'91| , [BK99| , |MS01| , [FS'1'02| ]. Such unitary 



transformations are obviously not unique, and since for most purposes it suffices to work 
with the projectors we hence refrain from using the unitary operators here. 

In view of the fact that is an orthogonal projector on the Hilbert space L^(R°') (8>C", 
one can ask if it is possible to satisfy also the relation (|2.5| ) exactly. In other words, to 
what extent can Vfj, be related to a spectral projection of 7i? (See ||HS88| , |CorOO| , |Cor01|| 
for examples.) We want to illustrate this question in the case where the principal symbol 
Hq of H possesses two well-separated eigenvalues < A^+i with constant multiplicities 
ki, and K+i, respectively, among the eigenvalues Ai, . . . , A/. For / = 2 this is exactly the 
situation that occurs in the case of a Dirac-Hamiltonian that we will discuss in some detail 
elsewhere pG|| . We also assume that there exists A G M separated from the spectrum 
spec(H) of H along with a fixed compact subset W C T*M'^ such that 



\ — \^{x,^) > Cm{x,^) and A,^+i(x, ^) — A > C"m(x, ^) 



(2.12) 



for all (x, ^) e T*M'^ \ W. It follows from these assumptions that for (x, ^) G T*W \ W one 
can replace the contour r< := lj)l=i in 

P<(x,0:=X^Pm(^,0 = ^ / Q{x,^,z)dz, (2.13) 

see ( ^.101 ), by a straight line r+ := {2; G C; z = \ + it, t G M} that avoids the eigenvalues of 
the principal symbol Hq as well as the spectrum of TC. Correspondingly, r> := lj|i=i/+i 
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is deformed into r_ given by r+ with reversed orientation. Thus, for (x,^) G T*M.'^ \ W 



If this relation held true for all (x, ^) G T*]R'^, = op^[P<] would be semiclassically 
equivalent to the spectral projection of H onto the interval (—00, A) given by 



whereas P> would correspond to 1(a,oo)('^)- For (x, ^) G W, however, it might happen 
that F-t crosses an eigenvalue of Hq. But the contribution to P<(x, ^) coming from the 
region W, where the eigenvalue functions have no sufficient distance from A, can be shown 
to be semiclassically small. Therefore: 

Proposition 2.2. // the eigenvalues Ai, . . . , A; of the principal symbol Hq are separated 
according to (HO) and the condition ( ^.ISj ) is fulfilled, the almost projection operators 



V< := op'^[P<], whose symbols are defined in ^2.1'^ , can be semiclassically identified with 
the spectral projections 1(_oo,a)('^) o-nd 1(a,oo)('^) of the operator Ti to the intervals (—00, A) 
and (A, 00), respectively. This means 

||P< - l(-oo,A)(H)|| = 0{K^) and \\V, - l(A,oo)(H)|| = 0{K^). 
A corresponding statement holds for the related orthogonal projectors V^, 

\\V, ~ l(-oo,A)(H)|| = and \\V, - l(A,oo)(H)|| = 

Proof. We use that fact that the contour of integration in the definition of P< in (p^.l3D 



can be deformed into T± outside the compact region W C T*R . To cut off the region W 
we choose a smooth and compactly supported function x ^ C^(T*M°') equal to one on W 
and use the corresponding partition of unity, 1 = x + (1 — x), to write 

(27ri)P<(a;,0~x(a:,0# / Q{x,^, z)dz + {1 - x{x,0)i^ [ Q{x,^,z)dz. 

In the second contribution, whose support is contained in T*]R'^ \ W, one can replace F^ 
by F-I-. We are thus left with the first term, which represents a symbol p{x,C,) in S°i(l) 
with compact support suppp C supp x- Here we apply a translation on T*]R'^ mapping 
p{x,^) to p{x,^) := p{x — Xq,^ — ^0)? with Xq and ^0 chosen such that ^ = is no longer 
contained in the support of p, i.e., p{x,0) = for all x G T*M'^. Using the fact that the 



Weyl operators corresponding to p and p are unitarily equivalent, see ||Fol89|| , we therefore 
obtain 

l|op^[p]|| = ||op^[p]||. 
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In order to estimate op^[p] in operator norm, we consider its action on a function ■?/' G 
y(R'^) (g) C" given by 

(op^[p]^)(x) = I j et(--^)^p(^,e)^(y) dyde. 

Since the symbol p vanishes in a neighbourhood of ^ = 0, one can perform an integration 
by parts after having inserted the operator (i^|^|~^(^ ■ dy))^ , N eN. We thus obtain that 
op^[p]'^ vanishes up to terms of order h°° and, since ^{M.'^) (8>C" is dense in L^(M'^) ® C", 
conclude 

\\op^[p]\\ = \\op^[p]\\=0{h°^). 
To finish the proof we now have to estimate 

/ (H-z)-' dz-^ [ op^[(l-x)#g(x,^,^)] dz 
in operator norm. To this end consider for 2; G r± 

((7^ _ _ op'^[i - x] op'^[Q(x, e, z)]){n -z) = {i- op^[i - x]) + 0{hn = 0{hn, 



which holds since x has compact support, and thus op'^fx] can be treated in the same 
manner as op'^[p] above. At this point the proof is complete, since for z G F-t the operator 
(7i — z) is invertible and its inverse has a norm that exceeds 0{h°°). □ 



3 Invariant algebra and Egorov theorem 

In this section our aim is to identify a suitable class of operators that is left invariant by the 
time evolution. Recalling the reasoning from the beginning of section |^, we are interested 
in a subalgebra of S^(l) whose filtration is respected by the time evolution generated 
by the one-parameter group U{t) = exp [—jT-Ct), where H is an essentially selfadjoint 
pseudodifferential operator with symbol H in the class S[!i(m). The following assumptions 
on the symbol H guarantee the essential selfadjointness of on S^iM.'^) ®C" (see |pS99||): 



(HI) H G S°i(m) is hermitian, 

(112) i^o + i is elliptic in the sense that \\{Hq{x,^) +i)~"'^||nxra < cm{x,C,)~^. 

Under the assumptions (HI) and (H2), U{t) therefore defines a strongly-continuous unitary 
one-parameter group. 

We now consider the time evolution of an operator B G B{L'^{M.'^) C") given by 

Bit) ■.= U{t)*BU{t), (3.1) 
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which is, of course, a bounded operator on L^(R ) ® C". According to Proposition IL^ the 



boundedness of B is guaranteed by choosing B G S^i(l). Moreover, a conjugation of (|3] 
with ^^=1 = idj;^2(Kd)0cn +0{h°°) results in a bounded operator so that 

/ I 
Bit) = J2 V^e^'^^^'Be^r.-HV.Ap^ = ^ e^^^"*P^i3P^e-^«^''* mod 0(n°°) (3.2) 

V,^=l U,fJ,= l 

in the operator norm. Here we have used the property e~T^'^^Vp = e'Ti'^^'^^'P^ modulo 
0{H°°) that follows from the Duhamel principle. Now, the principal symbol^ of TiV^ 
is a scalar multiple of the identity in the eigenspace P^^C'^ of Hq corresponding to A^, 
i.e., HqP^^o = X^Ppifi. Thus, ioi fi = u the operator exp (■^Ti'P^t) ;Bexp (^—j^HVut) is a 
pseudo differential operator with symbol in the class S''(l), see |[lvr98| , PGOOj . But when 



/i 7^ z/ the corresponding expressions are semiclassical Fourier integral operators. In that 
case the semiclassical limit at time t 7^ is different in nature from that at time zero. For 
a Dirac-Hamiltonian this phenomenon is related to the so-called "Zitterbewegung" which 
we will discuss in more detail in ||BG|| . Therefore, we are here interested in operators B 
with symbols in _B G S^[(l) for which U*{t)BlA{t) is again a semiclassical pseudodifferential 
operator with symbol B{t) G S^j(l). We hence introduce the following notion: 

Definition 3.1. A symbol B G S^[(l) is in the invariant subalgebra Sj^y(l) of the algebra 
S^(l), if and only if for all finite t the (bounded) operator B{t) = W{t)BU{t), B = op^[B], 
is a semiclassical pseudodifferential operator with symbol B(t) G S^i(l), i.e., 

S~ (1) := {B G Sl,il) ; symb^[ZY*(t)SW(t)] G S^i(l) for t G [0,T], g G Z} . 

This means that the invariant algebra S^^y(l) has a filtration, induced by the filtration 
of S^(l), which is invariant under conjugation of the corresponding operators with U{t). 
Due to the results of ||BGOO|| we expect that operators which are block-diagonal with respect 
to the projections are in the associated invariant operator algebra. This statement is 
made precise in Theorem |3.2| which is a variant of the Egorov theorem ||Ego69|| for general 
hyperbolic systems. 

Let us first consider an operator B with symbol B G S^[(l) that is block-diagonal with 
respect to the semiclassical projections, i.e., 

I 

According to the Heisenberg equation of motion ( p.2|) its time evolution B{t) is governed 
by 

^^B{t)r^^[H,B{t)]#. (3.3) 



•^We remark that before transferring operators to symbol level one can replace "P^ by Vf^ and employ 
the classical asymptotic expansion of the symbol . This will only amount to an error of order h°° . 
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Suppose now that B(t) has a (formal) asymptotic expansion 

oo 

j=0 

and use the composition formula of Lemma |1.2| together with the fact that the block- 
diagonal form of an operator B is preserved under the time evolution, see ( p.2| ). On the 
symbol level the diagonal blocks VyB{t)Vy then obey the following equation: 

j=0 

oo 

i,j=0 |a| + |/3|>0 



Here we introduced the notation F^g] := d?d^F for F e C°°(T*M'^) O M„(C), as well as 



7(a,/5) := 



(/3) •- 



2l-l+l/3||a|!|/3|!' 

oo 



J' 

oo 



j=0 

One hence has to solve, by taking [H^^q, B{t)i,y^t^ = into account, 

0<i<n-l 
j+\a\ + \l3\=n-l+l 

(3.4) 

Upon multiplying this commutator equation with the projection matrices P^^o from both 
sides one first realises that it is only solvable, if the diagonal blocks of the right-hand side, 
that we denote by Rn^u{t), vanish. The off-diagonal blocks on both sides of the relation 
(O) then yield the general structure of the solution, which reads 

B{t)uv,n+i = -P^t,o-B(t);yi,,„+iP^_o + Y^ — ^\ i'^-'^) 
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see also [ por95|| . This demonstrates that one obtains the off-diagonal parts of B{t)uy^n+i 
with respect to the projection matrices P^^o from the preceding coefficients of the asymp- 
totic expansion of B{t),,^. The diagonal parts then have to be determined by the condition 



that the commutator equation (|3.4|) must possess a (non-trivial) solution with initial value 
B{t)uu,n+i\f=o ~ -Bi/i/,n+i- Starting with n = 0, where the sum in ( |3.4|) is empty, one has to 
solve 



P, 



' d 1 



Pufi = 0. 



Expressions of this type have already been considered in ||SpoOO|| , where it was shown that 
the above equation is equivalent to 



d . ~ 

see also appendix Here we have defined the hermitian n x n matrix 



(3.6) 



(3.7) 



according to ( |A.4| ) and (|A.5| ) of appendix 0. Now, equation ( p.6|) is trivially fulfilled for 
z/ 7^ /i, and the case u = fi has already been considered in |p.vr98| , |BN99| , |BGOCI|| , where it 
was shown that the solution reads 

B{t)^^^o{^,x) = d-l{x,^,t)B^^^o{<!>l{x,0)duu{x,^,t). 

In this expression $^ : T*M.'^ — > T*M.'^ denotes the Hamiltonian flow generated by the 
eigenvalue Aj, of Hq, and the transport matrix di^i, is determined by the equation 



duu{x,^,t) + iH^i,^i[^l{x,^))duu{x,^,t) = 0, dyi,{x,^,0) = idi 



{3.1 



One has thus fixed the coefficients Bit)^^^ = P^fiB{t)QP^^o, i.e., the principal symbol of 
B{t), since the off-diagonal terms B{t)u^fl = PyfiB{t)QP^fi vanish and therefore trivially ful- 
fill (|3.4|) . According to ( p.5|) we hence have also determined the off-diagonal parts of the sub- 
principal term B{t)^^^i, which vanish as well. The diagonal contributions P^flB{t)yuiP^Q 
with respect to the projection matrices obey [P,,,0) P^lfiB{t)yy^lP^^Q\ = and thus can be de- 
termined from the relation ( ^.4[ ). As in ||ivr9(^ , |BGOO[| , we hence obtain a recursive Cauchy 
problem for the coefficients B{t)uu,n and are now in a position to state: 

Theorem 3.2. Let H G S°i(m) he hermitian with the property 



^.■S(a;,Ollnxn<Q^ for all {x^Oen" and\a\ + \(3\^j>2-5_ 



'jO, 



(3.9) 
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and such that the conditions (HO) and (H2) are fulfilled. Furthermore, suppose that B G 
S^j(l) is block- diagonal with respect to the semiclassical projections defined in i\2.!^ ), 

I 

Then B is in the invariant algebra S?^^{1) introduced in Definition \3. 1\ , i.e., B{t) is again 
a semiclassical pseudodifferential operator with symbol B{t) G S^i(l). Furthermore, the 
principal symbol of B{t) is given by 

I 

B{t)o{x,0 = J2dUx,C,t)B,,4<!>i{x,0)duAx,^,t), (3.10) 
1^=1 

where is the Hamiltonian flow generated by the eigenvalue \y of Hq, andd^i, is a unitary 
n X n matrix which is determined by the transport equation l \3.^ . 

Proof. As in | Ivr98 , BGOO | we start by rewriting ( p.4|) for the diagonal block of B(t)^i,^n 
with respect to P^^q in the form 

^ [rf;^^(x,e, -t)(P^,o5(t)..,„P^,o) o <i>;*^-(x, 0^.^(3;, -t)] 

0</<n-l 
j+\a\ + \l3\=n-l+l 

(3.11) 

where ci^^ is determined by the transport equation 

d,^{x,i,t)+iH,^,^i{^'^'^-{x,i))d,^{x,i,t), d,^(x,e,0) =idcn, 

that generalises ( p.8|) also to the off-diagonal transport. And since H^^ i is hermitian, the 
solution dvfi is a unitary n x n matrix, which in the case u n is obviously given by 

d,^(x,e,t) = e-^^-^(^'«)*. 

In order to obtain estimates on the derivatives of the symbols Pfi^B {t) yy^nif) Pfifi one has 
to control the behaviour of the flow $^ generated by the eigenvalue Xi, of Hq. To this end 
we first notice that Hq G S(m) implies the bound |A,^(x, ^)| < cm(x, ^) on its eigenval- 
ues. Furthermore, due to the hyperbolicity condition (HO) the projections Py^ onto the 
eigenspaces of Hq are in S(l). We then consider the first order derivatives (|a| + \(3\ = 1) 
of the relation 

Hq{x,OPu,o{^,0 = A.(x,OP.,o(a:,0, 
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which exist since the eigenvalues are smooth functions on the phase space T*M'^, see 
equation (|2.8| ). One thus obtains 



Now, since Pu,o{x,OPu,o ^^]{^^OPi^,o{x,0 = 0, a multiphcation of the above equation with 
Pu,o{x,C,) from both sides yields 



and hence 



A,g(x,OP.,o(2;,0 = P.,o(a;,0^o(a;,OS^.,o(x,0, (3-12) 



|\ (")p 11 —A\p (")p II ^;;|lf7'(")|l 



Hq G S(m) therefore implies that the first order derivatives of Xi, are bounded by the 
order function m. One can continue this argument by successively differentiating equation 
( p.l2| ), and concludes that X^ idc" G S(m) for all z/ = 1, . . . ,1. In particular, the property 

||^oS;;(x,0|Lx„<C.,^ for |a| + |/5|>l, 

which follows from (|3.9|) , transfers to a corresponding growth property of the eigenvalues 
of ifo: 

|A.g(x,0| <C,,/3 for \a\ + \P\>l. 

Therefore, the Hamiltonian flows $* exist globally on T*M'^ such that < C^,/? 

for all a, /? G ^^id for all finite times t G [0, T], see [|Rob87|] . This property guarantees 
that i? o <|)^ G S(l) for all B G S(l). Concerning the unitary matrices d^^^ the following is 
true: 

Lemma 3.3. If the subprincipal symbol Hi ofTi satisfies \\Hi |^j||nxn ^ Co,/3 for all \a\ + 
\(3\ > 1, then Wd^,^ P){x,^,t)\\nxn < C'^^p for allt G [0,T], |a| + > 1 and u, = 1, . . . ,1. 

For the proof of this lemma see [[lvr98|| . With these properties at hand one can integrate 
equation and solve for P^fiB{t)yi,^nP^ifi by conjugating with d^^{x, ^, —t) and shifting 

the arguments by ^i"""^ (which only amounts to an actual shift in the case v = jj). For the 
principal symbol of B{t) one thus obtains 

5(t),,,o(x,0 = d,,{^'^{x,i),-t)B,,^o{^'^{x,i))d-l{^l{x,i),-t), 

which is the only block of B{t)uu,o with respect to P^^, = 1, . . . , /, that is different from 
zero. Using 

= <(a:,e,t) =<^(x,e,t), (3.13) 
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see [ |BN99|| , one finally obtains ( p.lO| ). For the higher coefficients B{t)yy^n, n > 1, one 
employs the Duhamel principle and uses that fact that the sum in is taken over 

indices with \a\ + \13\ + j > 2, and thus involves terms in S(l), in order to conclude that 
B{t)yy^n £ S(l). This shows that one has found an asymptotic expansion in S^i(l) for the 
symbol of lA*{t)BU{t) that can be summed with the Borel method to yield a complete 
symbol. □ 

This theorem shows that, for finite times t, one can associate to a (semiclassically) block- 
diagonal symbol B G S^i(l) a symbol B{t) G S^i(l) whose quantisation op'^[i?(t)] is semi- 
classically close to B{t) = U*(t)BU(t), i.e., 

\\B{t) - op^[B{t)]\\ = 0{K^) for all t G [0,T]. 



This is a semiclassical version of the Egorov theorem ||Ego69|| , which was originally for- 



mulated for the case of scalar symbols. A weaker version that is also sometimes referred 
to as an Egorov theorem (see, e.g., |PST02|| ) would only assert that one can evolve the 



principal symbol Bq of B into a symbol B{t)o, as given in ( |3.1CI| ), such that its quantisation 
op^[B(t)o] is ^-close to the time-evolved operator B(t), i.e.. 



||i3(t)-op^[i?(t)o] II =0(/l). 



This statement is clearly covered by Theorem p.2| , since the quantisation of the difference 
B{t) — B{t)o G S^j~^(l) yields a bounded operator with norm of order H, see Proposition 



We will now show (generalising results of Cordes ||Cor834 porOOl , |Cor01| ]) that the 



semiclassical block-diagonal operators exhaust all operators with symbols in the invariant 
algebra (1). 

Proposition 3.4. The invariant subalgebra of S^(l) consists of precisely those 

B G S^[(l) which are semiclassically block- diagonal with respect to the projections P^, fi = 
1, . . . , /, defined in l \2. 1 1\ ) of Proposition |^. i| , i.e., 



I 

i? G s?^, c s,7(i) ^ fi^^p^#i?#p^. 

Proof. Consider an operator B with symbol B G SJJ(I), whose equation of motion is given 
by ( p.3|) . For the symbol of the time-evolved operator we now assume an asymptotic 
expansion 

oo 

B{t) ~ ^-'^'mj 

in S^i(l). Furthermore, one can use ( p.7|) to separate ( |3.3| ) into blocks with respect to P^, 
fi = 1, . . . ,1. For the off-diagonal blocks {u ^ fi) one therefore obtains 

|p(t).,~-^[i/,P(t).,]#, (3.14) 
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where B{t)y^ := Pui^B{t)i^P^ ~ Yl'jLo^ '^^^ B{t)y^^j. In leading semiclassical order the 
factor hr^ on the right-hand side of equation (|3.14| ) enforces the condition 

[Hf), -B(t)i//i,o] = — ^ti)B{t)i,^^Q = 0. 
Since 7^ \u for ^ this immediately yields Bit),^^^ = 0. Furthermore, 



# 



Again the leading order on the right-hand side has to vanish, i.e., 

[Ho, = 0. 

This means that the symbol B{t)y^^i must be block-diagonal with respect to the projection 
matrices P^^ G S(l). But 

PA,o5(t).^,iPA,o = symb]r [h-^ {PxHB{t),,, - B{t),^,o)#Px)] = 0, 

since Bit)^^^ = for z/ 7^ /i. Iterating the above procedure we see that if i? G Sinv(l), then 
it has to be block-diagonal with respect to P^, /i = 1, . . . , /. This proves one direction 
asserted in the proposition. The other direction, that the block-diagonal operators form a 
subset of the invariant algebra, is contained in the Egorov theorem ^.2|. □ 



At this point we want to add a comment on the transport equation ( p.8| ) that not only 
occurs in connection with an Egorov theorem, but also in a WKB-type framework. In this 
context Littlejohn and Flynn [|LF91|| introduced a splitting of the analogue to Hy^ i (defined 
in equation (|3.7|) ) into two contributions, one of which is related to a Berry connection 
|Ber84|| . Subsequently Emmrich and Weinstein ||EW96|| generalised the approach of |P^F91 



and gave a geometrical interpretation for the second contribution, which they related to a 
Poisson curvature. We now want to identify the two contributions in the present situation, 
i.e., in To this end we calculate Hy^i = Pu,iHo + Pu,qHi + ^{Py^, Hq} using 

—PuflPu,lPufi + (1 — Pufl)Pv,l{)- — Pufl) = -^i^ufi, Pu,o}, 



which follows from the condition P^jj^Py ~ Py and the composition formula in Lemma 
Thus 



A i ' 

PvfiHylPyQ = PyfiHlPyQ + i — Pyfi{Py^O, Pyfi}Py^O ^ " \Pu,o{Pu,0, P'n,o}Bu,0- 

ri=l 



The relation Pyfl{Pyfl, Prj,o}Pufi = —Pu,o{Prifi, Pr],o}Pu,o and the spectral representation 
Hq = A^P^^o now allow to rewrite the expression (|3.7|) for Hyy i as 

Hyy^l = Hy^ 

Berry 

+ Hy 

,Poisson 

+ PyfiHlPyfl 
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with 

,Poisson • „ ( ^uPu,o{Pufi-, Pufi}Pu,Q + Pufi{Pufii Hq — Xi,Pufi}Pu,i 

This corresponds exactly to the sphtting discussed in [[EW96|] , see also [^poOO|| . 



4 Dynamics in the eigenspaces 



According to the Egorov theorem |3.2| , the semiclassical calculus outlined above results not 
only in a transport of the principal symbols of observables by the Hamiltonian flows , but 
also in a conjugation by the (unitary n x n) transport matrices d,^iy The latter define the 
dynamics of those degrees of freedom that on the quantum mechanical level are described 
by the factor C" of the total Hilbert space L^(R°') ® C". Our intention in this section now 
is to develop combined classical dynamics of both types of degrees of freedom, i.e., those 
described by the Hamiltonian flows and those that are represented by the conjugations. In 
this context the fact that the conjugations enter along integral curves of the Hamiltonian 
flows introduces a hierarchy among the two types of degrees of freedom. 

In a first step we confirm that the dynamics represented by the transport matrices d^u 
take place in the eigenspaces of the principal symbol Hq in C". To this end we notice that 
since at every point (x, ^) e T*R'^ the projection matrices Pufi{x,^) yield an orthogonal 
splitting of C" and have constant rank k^, they define /cj^-dimensional subbundles '■ E'^ ^ 
T*W^ of the trivial vector bundle T*R°' x over phase space. The fibre E^^s,) ~ ^iJ^l^^O 
over (x, ^) G T*R'^ is given by the range of the projection, i.e., E'^^^-^ = Pi,^o{x,^)C"'. 
Furthermore, the canonical hermitian structure of C" induces a hermitian structure on the 
fibres. We now intend to interpret the conjugation by d^,^ as a dynamics in the eigenvector 
bundle and for this purpose notice: 

Lemma 4.1. The restricted transport matrices d^^{x,^,t)P,y^o{x,C) Provide unitary maps 
between the fibres E^^ and E^t ■ 

Proof. In order to see that di,i,{x,^,t)P^^o{x,$,) maps E'j^^^s^ into E^tf^^^^ we show 

P,,oiKix,0)du.ix,^,t)P,,oix,0 = rf..(x,e,t)P.,o(x,0- (4.1) 

This relation is certainly true for t = where both sides yield P^^q. Moreover, the derivative 
with respect to t of the left-hand side reads 

{A,, P.,o}($t(a;, 0)duu{x, t)P,,o{x, - ^P.fiWx, 0)H.u,i{K{x, 0)d.u{x, t)P,,o{x, 0, 



which equals 

-iH^uAKix, 0)Pufi{'^i{x, 0)duu{x, ^, t)P^,o{x, 0, 
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since the commutator [Pu,o, Hui^^i] can be calculated as (see equation (p^)) 

here we have used ( |A.2| ) and P^q = Py^. Thus, Pufl{^*y{x,^))duv{x^^,t)Pufl{x,S,) fulfills 
the same differential equation with respect to t as duu{x,^,t)P„^o, and this finally implies 



the validity of equation { ^TL 



In order to see the unitarity, one has to show that d,yi^{x, ^,t)P,y^Q(x, ^) is an isometry 
whose range is the complete fibre E^t (j, The first point is clear since d^i, is unitary on C" 
and the fibres inherit their hermitian structures from C". The second point follows from 
the observation that the transport provided by d^u can be reversed: Given v{^l{x,^)) e 
the vector P^^o{x,^)d,^i^{^^j^{x,^), —t)v{^l{x,^)) lies in E'^^.^) mapped to 

v{ii{x, 0) by d,,{x, e, t)P.,o(x, 0, see §J^. ' □ 

According to the above, the action of duu{x,^,t) on a section in E'^ can be viewed as a 
parallel transport along the integral curves of the fiow <l>^. If one now introduces sections 
of E" that yield orthonormal bases {ei(x, ^), . . . ,ek^{x,^)} of the fibres E^^^)^ repre- 
sentations of di,y{x,^,t) in these bases are unitary ky x ky matrices Dy{x,^,t). Since the 
principal symbol Hq of the operator 7i is hermitian (on C"), a preferred choice for the sec- 
tions {ci, . . . , Cfc,^} would consist of orthonormal eigenvectors of Hq. However, this choice 
is obviously not unique because it amounts to fixing an isometry Vy{x,^) : C^" ^(x^)^ 
such that Vy{x,^)V*{x,^) = Py^o{x,^) and V*{x,^)Vy{x,^) = id^k^. Here one still has a 
freedom to change the isometry by an arbitrary unitary automorphism of C'^". Having 
chosen an isometry Vy{x,C,) for every fibre E^^^^ in a smooth way, the n x n transport 
matrices dyy{x,^,t) are mapped to unitary k^ x k^ matrices 

D,(x,e,t) := V:{^iix,0)dUx,^,t)K{x,0. (4.2) 
Their dynamics follows from the transport equation ( ^.8| ) as 

Dy{x,^,t) + iHy{^i{x,O)Duix,^,t) = with D,(x,e,0) = idcfc., (4.3) 
where the hermitian k^ x k^ matrix is derived from ( |3.7| ) for /i = i^. 



Hu = -iy K*{^.,o, Pu,o}Vy + i{A„ v:}Vy + 

What is of more importance for later purposes than the non-uniqueness of this represen- 
tation, however, is the fact that the above construction allows to introduce a skew-product 
fiow over the Hamiltonian fiow thus refiecting the hierarchy of the two types of degrees 



of freedom. See ||CI''S82| ] for a definition of skew-product fiows and cf. | [BK99b| | where these 



occur in the context of a semiclassical trace formula for matrix valued operators. At this 
stage now provisionally consider 

: T*R°' X \J{ky) T*R"' x \J{ky), 
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defined by Y^{x,C,,g) '■= {^l{x,^), D,y{x,C,,t)g), which yields a flow on the product space 
T*M'^x U(A;^) due to the cocycle relation D^{x,^,t + t') = D^{^l{x,^),t')Dj{x,^,t). Later 
we are interested in ergodic properties of such skew-product flows, and these are indepen- 
dent of the particular choice of the sections {ei, . . . ,ek^}. Here we remark that in some 
cases the point of view advertised above might turn out too general. It can indeed happen 
that the flbre part of the skew-product flow does not require the complete group U{ki^). 
E.g., in [|BGK01|1 a situation was considered where k^, = 2j + 1, j G hN, and the trans- 



port matrices 0,^ were operators in a 2j + 1-dimensional unitary irreducible representation 
of SU(2). This fact could be identifled by the observation that when (x, ^) ranges over 
T*M'^, the skew-hermitian matrices iHy{x,^) generate a Lie subalgebra of u(2j + 1) which 
is isomorphic to su(2). 

In the general case one therefore should not necessarily expect that the transport ma- 
trices generate all of \]{ky), but only a certain Lie subgroup. In order to identify this 
group we consider the Lie subalgebra 

(i#,(x,0; (x,OgT*m^>cu(M (4.4) 

generated by the skew-hermitian matrices iHy{x,^). Via exponentiation of this subalgebra 
one hence obtains a Lie subgroup G C U{ky) that is compact and connected. To be more 
precise, the result of the exponentiation is a /cj,-dimensional unitary representation p of 
G. Its Lie algebra g then is embedded in ([4.4| ) via the derived representation dp. In this 
setting the transport matrices D^, are operators in the representation p, i.e., Dy{x,^,t) = 
p{gu{x,C,,t)). Hence we are now in a position to deflne the skew-product flows 

: T*R'^ xG^ T*R'^ x G (4.5) 

through 

Y^{x,^,g) = {<l>i{x,0,gu{x,U)g)- (4.6) 

These flows leave the product measure da;d.^d5f on T*M.'^ x G invariant, which consists of 
Lebesgue measure dx on T*M'^ and the normalised Haar measure dg on G. Moreover, if 
one restricts the Hamiltonian flows to compact level surfaces of the eigenvalue functions 
at non-critical values E, 

n,,E := K\E) = {(x,0 G n"; K{x,0 = E}, 

the restrictions of the skew-product flows to Qi^^e x G leave the measures di{x,^)dg 
invariant, where di{x,^) denotes the normalised Liouville measure on Q^^e- 

Below we are interested in the question under which conditions imposed on suitable 
classical dynamics quantum ergodicity holds, see section ||. In analogy to ||BGOO|| one 



approach to this problem would be to consider the restriction of the skew-product flow Yj; 
to Qiy^E X U(fci/): its ergodicity with respect to the product measure that consists of Liouville 
measure on ^ and Haar measure on U(fcj,) implies quantum ergodicity. Since, however, 
the dynamics in the eigenspaces is completely flxed by a restriction to the group G, the 
dynamical behaviour of the flow Y^ is determined by that of Y^. One hence concludes that 
in order to proof quantum ergodicity one requires the following condition (see Remark |6.3|): 
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(Irr,^) The representation p : G ^ ^{k^) is irreducible. 

In the sequel we always assume this to be the case. 

Our intention now is to relate the dynamics in the eigenspaces, given by the conjugation 
with the transport matrices (1^,^, to proper classical dynamics. To this end we require a 
symplectic manifold with the dynamics realised in a Hamiltonian fashion. For this purpose 
Kirillov's orbit method ||Kir76|| provides the necessary tools: it relates the unitary irre- 
ducible representation (p, C'^") to a coadjoint orbit O of G, which is a symplectic manifold. 
Moreover, the conjugation dynamics is realised in terms of the coadjoint action of G on 



O. As in the case of G = SU(2) considered in [ [BGK01|| , this setting then also allows to 



introduce a Moyal-type quantisation such that hermitian matrix valued symbols can be 
uniquely related to real valued functions on the symplectic product phase space T*M'^ x O. 
Let us now recall some properties of coadjoint orbits ||Kir76 |: The adjoint representation 



Ad : G aut(g), g t-^ (Tel)((yf), of a Lie group G on its Lie algebra g = TgC is defined 
as the differential Tgl of the inner automorphism l{g) : G —>■ G, x ^ gxg~^, g E G, at 
the identity e E G. The coadjoint representation of G on the dual Lie algebra g* is then 
provided by the dual Ad* := (Ad^-i)* of the linear map Ad^-i, i.e., 

(Ad;(A),X) = (A,Ad,-iX), 

for X E Q and A G g*; here ( , ) : g* x g ^ M denotes the dual pairing between the vector 
spaces g and g*. A coadjoint orbit Ox through A G g* then is an orbit of this group action, 

O, := {Ad;(A); geG}c g*. 

If G is compact, 0\ is a smooth embedded and compact submanifold of g*. One of the 
main features of coadjoint orbits is their symplectic structure |Kir76|. 



Proposition 4.2. Let G be a connected Lie group and O d g* a coadjoint orbit. Then O 
is a symplectic manifold and there exist unique symplectic forms on O such that 

a±(A)(ad;^A,adyA) =±(A, 

for all X E O and X,Y G g. Here ad* denotes the differential of the coadjoint action and 
[ , ] is the Lie bracket on g. The forms cr^ are referred to as the coadjoint orbit symplectic 
structures. 

Furthermore, let Gx '■= {g G G; Ad* A = A} denote the isotropy subgroup of A G g* 
under the coadjoint action. Then this is a closed subgroup of G, and so the quotient G/Gx 
is a smooth manifold with smooth projection tt : G —>■ G/Gx such that one can identify 
G /Gx = Ox via the diffeomorphism k : gGx ^ Ad* A. Moreover, since the coadjoint action 
on Ox preserves its symplectic structure and is obviously transitive. Ox is a symplectic 
homogeneous space. In the following we only need one symplectic structure that turns Ox 
into a symplectic homogeneous space and therefore now fix cr := cr+. 

Our next goal is to construct a certain quantisation of the symplectic manifold Ox- On 
the classical side one considers suitable functions on the phase space Ox as observables; 
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here we choose functions that are integrable with respect to the volume form drj that 
arises as the maximal exterior power of the symplectic form a. A Hamiltonian dynamics 
is then generated by a smooth real valued function h on 0\ through the association of a 
Hamiltonian vector field Xh according to (T{Xh, ■) = dh. On the quantum side observables 
are hermitian endomorphisms of the representation space V with inner product (■, ■)y. 
A Moyal quantiser now assigns to a hermitian A G C{V) a function a e L^{0\) such 
that p{g)Ap{g~^) is mapped to a o Ad*-i. This covariance property then ensures that the 
dynamics given by the conjugation with D^, = p{gu) is represented on the phase space Ox 
through the coadjoint action of g,y. Quantisations of this type were constructed by Simon 
Sim80|| , who introduced suitable Berezin symbols representing A G jC.{V). However, here 



we will closely follow ||FGV90|| , where it is demonstrated that one can obtain a quantisation 
with an additional tracial property that turns out to be useful later. 

In the present context G is a matrix Lie group, i.e., a closed subgroup of GL.„(C), and 
its Lie algebra g is a subalgebra of M„(C) with the matrix commutator as Lie bracket. Thus 
there also exists a non-degenerate symmetric bilinear form B on g, given by B{X,Y) = 
Re tr(Xy), that allows to identify g with its dual g*. Consider now the unitary irreducible 
representation (p, V) of G and fix a highest (real) weight A G t* corresponding to this 
representation, where t is the Lie algebra of a suitable maximal torus T G G. Since one 
can identify g and g* via the bilinear form B, one can regard the highest weight as A G g*. 
Up to a phase, to this highest weight there corresponds a unique normalised weight vector 
wx G V. Now define the map J : V ^ g* hj 

{Av),X) := {v,dp{X)v)v, 

such that in particular J{wx) = A. This map is equivariant in the sense that J{p{g)v) = 
Ad* J{v). Since the weight space of the maximal weight is one-dimensional one obtains 

j-i(A) = {zwx; z eC, \z\ = 1}, 

and therefore 

J~\Ox) = {pig)wx; geG}. 

This setting now allows to associate to points t] E Ox vectors f ^ G V that are unique up 
to a phase. For this purpose one chooses a measurable section t] g^, in G ^ G/Gx with 
Ad*^(A) = 7] and gx = e, which is possible due to the fact that Ox = G/Gx is an orbit of 
the coadjoint action. Then define for every r] E Ox a. vector t>^ := p{gri)wx, which can also 



be viewed as a coherent state |Per86|| . The equivariance of the map J now implies that 
J{v^) = 1], such that is unique up to a phase. This finally allows to define for every 
A G C(y) the unique covariant symbol 

Qa{v) ■■= {v^,Av^)v- (4.7) 

In fact, Qa '■ Ox ^ C is continuous. We denote the space of covariant symbols that are 
constructed according to the above scheme by 5a := {Qa', A G C(y)} and recall from 

IFuvg 
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Lemma 4.3. The map C{V) Sx defined in equation i\4- Tj ) is one-to-one. 

Now consider a normalised vector w & V and the associated orthogonal projector 11^ 
onto Cw C V. Since, in the language of quantum mechanics, expectations of an observable 
A G C{V) in the state w read {w,Aw)v = tr(y4n„), one would like a Moyal quantisation 
to represent tT{AB) as 

QA{v)QBiv) dr/. 



o 



A 



This relation, however, does not hold. Considering C{V) as a (finite-dimensional) Hilbert 
space with inner product tY{A*B), we are hence looking for an isometry CiV) L'^{0\). 
To this end one notices that S\ being a finite dimensional subspace of L?{Ox), the Riesz 
representation theorem ensures for every A G CiV) that the linear form : S\ ^ £- given 
by La{Qb) '■= tT{A*B) can be represented in terms of a unique function P4 G Sx such that 

La{Qb) = tT{A*B) = [ Pa{v)Qb{v) dr]. 



A 



Since according to Lemma ^73| the spaces jC{V) and Sx have the same (finite) dimension, the 
map C{V) 3 A Pa E Sx is a (linear) bijection that can as well serve as a symbol map; 
Pa is then called contravariant symbol. We remark that in order to satisfy the natural 
condition Pjd^ = 1, the volume form on Ox has to be normalised such that vol{Ox) = 
dimV. Both the covariant and the contravariant symbol of A G C{V) fulfill the covariance 
condition 

Qpig)Ap{g-^){v) = <5A(Ad*-i 7]), 
Pp(9)Ap{g-')iv) = ^A(Ad;„i 7]), 

for all g & G and all r] G Ox- However, in order to obtain the desired isometry from C{V) 
into L'^{Ox), one is forced to introduce a symbol map that in a certain sense lies in between 
Q and P. 

In |[FGV9CI|| it is shown that the operator K on Sx that maps Qa to Pa is bijective and 



positive. It therefore allows for a (positive) square-root K^^'^ which can be used to define 
a symbol map with all desired properties: 

Definition 4.4. For A G C{V) the Stratonovich-Weyl symbol symb'^^[y4] G Sx is given 
by 

symb^^[A] := K'^'Qa = K'^'^Pa- 
Summarising the above finally yields [|FGV90|| : 



Proposition 4.5. The symbol map A ^ symb'^^[A] has the following properties: 
(i) It is a linear one-to-one map from C{y) to Sx, 
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(ii) symb^^[A*] = symb^^[A], 
(ill) symb^^[idy] = 1, 

(zv) SYmh''^[pig)Apig-')]{7^) = symh'"^ [A]{AdU v) for all r^eOx, geG, 



(v) [ symb'^^[A]{ri)symb'^^[B]{ri) dri = ti{AB). 

J Ox 



In order to make the relation between A G C{V) and its symbol symb'^^fA] explicit, 
one introduces a (hermitian) quantiser : Ox '^(V) such that 

symb^^[A] = tr(AAA) and A=[ symb^^[A](r/)AA(77) d?]. (4.8) 

J Ox 

As shown in ||FGV90|| , the quantiser can be expressed in terms of generalised spherical 



harmonics associated with those unitary irreducible representations that appear in the 
regular representation {j{g)f){vi) = /(Ad*_i rj) of G on 5a. 

With this formalism at hand one can now transfer the dynamics of a (hermitian) B G 
C{V) given by a conjugation with D{t) = p{g{t)), B ^ B{t) = D-\t)BD{t), to the 
coadjoint action of g{t) on the symplectic manifold Ox via the relation symb'^^[i?(t)](?7) = 
symb"^^ [5] (Ad*(j-) rj). The symplectic structure on Ox defined by the form a, furthermore, 
allows to identify the dynamics r] i— > Ad*^-^-) rj as being Hamiltonian. To see this assume 
that D(t) is determined by 

D{t) + iHD{t) = with D{0) = idy, (4.9) 

where H G C{V) is hermitian; compare equation ( [4.3|) . On the one hand now, a Hamilto- 
nian flow 7] \^ ri{t) can be introduced on Ox that is generated by the Stratonovich-Weyl 
symbol of H. The associated Hamiltonian vector field '^symb^^ [h] is then defined through 

^i.KjmhS^[H]^ ■) = dsymb^^[i/], 

so that the time evolution f(t){ri) = f{rj{t)) of a function / G C°°(Ca) is governed by the 
equation 

fit) = {symb^^[ff],/(t)} = = -df{t){X,^^^sw^^^). (4.10) 

On the other hand, differentiating Ad*(-^-) rj with respect to t yields 

^Ad*(,)r^|^^^ = -ad;,^r^, (4.11) 

where Xh G g is the generator of the curve g{t) in G which, according to equation ([4.9| ), 
is related to H via dpiXn) = —iH. A comparison of ( ^.lOj ) and ( [4. 11] ) then shows that 
the dynamics provided by the coadjoint action Ad*(-() coincides with the Hamiltonian flow 
generated by the symbol symb'^^[if]. 
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As an ultimate outcome of the above formalism we are now in a position to introduce 
a skew-product flow on the symplectic phase space T*M'^ x 0\ that completely determines 
the time evolution of the u-th diagonal block of an observable on the level of its principal 
symbol. Explicitly, this flow is given by 

: T*R'^ xOx^ T^M'^ x Ox (4.12) 

with 

:= {K{x,0,Mli.,^,^v): (4-13) 

it leaves the product measure dxd^ drj invariant. 

Consider now a semiclassical pseudodifferential operator B with symbol B G S^(l). 
Mod 0{h°°) the quantum dynamics preserves the diagonal structure of its blocks VyBVu. 
According to the Egorov theorem |3.2| , together with the definition ( [1.2| ), the principal 
symbol of VyB{t)Vu hence reads 

od:{x, e, t) (KT^oK) o)D.{^. e, t)v:{x, o. (4.14) 



We now exploit the possibility, explicitly provided by (|4.^ ), to uniquely represent the value 
of K*^oK : T*M'^ C{C^") in terms of a Stratonovich-Weyl symbol, 

6o,.(x,e,r/) := symb^^[(V;*5oK)(x,0](^)- (4.15) 

The dynamics of the principal symbol in this representation is now summarised in the 
following variant of the Egorov theorem: 

Proposition 4.6. The Stratonovich- Weyl symbol b{t)o,u associated with the principal sym- 
bol of the operator VyB{t)Vy is the time evolution of bo^^ under the skew-product flow 
defined in equations ^.1^} - ^.!^} , i.e., 

6(t)o,^(x,e,r?) = bQ,^{Yl,{x,i,T])). 

Proof. According to (^4.14|) and ( [4. 151) , b{t)o^i, is given by 

6(t)o,.(x, V) = symb^^ [pig-\x, e, t)) {V:BoV,) (<l>t(x, 0)pi9A^. ^, t))] (r^). 



which due to the covariance property (iv) of Proposition [4.5| reads 

6(t)o,.(x,e,r7) = symb^^[(K*5oK)($t(a:,0)](Ad;(,,^,)r/) 
= 6(t)o,.($t(x,0,Ad*^(,_^,,) r]). 



□ 
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5 Trace asymptotics and a limit formula for averaged 
expectation values 

A fundamental ingredient in the asymptotics of eigenvectors we are aiming at is a semi- 
classical limit formula for the expectation values of bounded operators B on L^(]R'^) ® C*^. 
Below we will obtain a Szego-type formula which connects semiclassically averaged ex- 
pectation values with objects that can be calculated from the principal symbol Bq of the 
operator B and therefore allow for a classical interpretation. On the so defined classical 
side we fix a value E for all eigenvalue functions X,y, z/ = 1, of the principal symbol 
Hq with the following properties: 

(H3,y) There exists some e > such that all X^^{[E — e,E + e]) C T*]R^ are compact. 

(H4,^) The functions X^, shall possess no critical values in [E — e, E + e]. 

(H5,^) Among the level surfaces Q,^^e = u = 1, ... ,1, at least one is non-empty. 

In addition to (HI) and (H2), which imply the essential selfadjointness of the operator 



Ti, these conditions guarantee as in the scalar case |PS99|| that for sufficiently small h the 



spectrum of Ti. is discrete in any open interval contained in [E — e,E + e]. This setting 
now allows us to generalise the constructions made in ||BGOO|| to Hamiltonians with non- 
scalar principal symbols: The expectation values of an operator B will be considered in 
normalised eigenvectors ipj of Ti with corresponding eigenvalues Ej in an interval I{E, K) = 
[E — hw, E + hw], LJ > 0, such that I{E, K) C [E — e, E + e] ii h is small enough. On the 
classical side the Hamiltonian flows $^ generated by the eigenvalue functions X^ will enter 
on the level surfaces VL^^e- Regarding these we assume: 

(H6i/) The periodic points of with non-trivial periods form a set of Liouville measure 
zero in VLu^e- 

The quantities appearing on the classical side of the limit formula turn out to be averages 
of smooth matrix valued functions B G C°°(T*M'^) ® M„(C) over VtyE with respect to 
Liouville measure, for which we introduce the notation 



The main result of this section is now summarised in the following Szego-type limit formula: 

Proposition 5.1. Let Ti he a semiclassical pseudodifferential operator with symbol H G 
S°i(m), such that the principal symbol Hq satisfies the assumptions (H0)-(H2) and (H3y)- 
(H6y) for all u = 1, . . . , / . Furthermore, let B be an operator with symbol B G S°j(l) and 
principal symbol Bq. Then the limit formula 

holds. 
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Proof. Adapted to the spectral localisation mentioned above we choose a smooth and 
compactly supported function g G C^(M) such that g{\) = A on a neighbourhood of 
[E — e, E + e]. Furthermore, we apply the semiclassical splitting of the Hilbert space 
L'^iW^) C" given by the projection operators Vu, 



C" = ran Pi 



ranVi mod h° 



(5.2) 



and the corresponding decomposition Ti. = Ylu=i'^'^'^ (mod 0{h°°)) of the Hamiltonian. 
By employing the generalisation of the Helffer-Sjostrand formula to matrix valued operators 
developed in |pim93| , pim98|| , we represent giTi.) = Yl^y=ig{T^'Py)Vu (mod 0{h°°)) with 



g{nv,)v. 



d,~g{z){n 



n d^, 



where g is an almost-analytic extension of g. Since the principal symbol HqP^ q of TiVi, is 
scalar, HoPi^fi = XuPu,o, when considered to act on sections in the eigenvector bundle E'^, 
one can use the methods of [pS99|| to show that on X~^{[E — e,E + e]) the asymptotic 
expansions of symh^ [gCHVu)] and of symh^ [HV^] coincide. Below we will always employ 
the spectral localisation to the interval I{E,h)^ and since SYmh^[g{HV,y)] G S°i(l), one 
can therefore now assume that H G S°i(l). Furthermore, the decomposition ( [5. 21 ) allows us 
to employ the techniques of |PS99|| in the same manner as in PGOO|| . Hence, if x ^ 
with X = 1 on I{E, K) and supp x C [-E — e, -E + e], the operator 



I 

X(H)5^P. = $^e-^^^''*x(^n)n mod 0{h° 



u=l 



has a pure point spectrum. Moreover, each of the operators e^T^^^'^^xO^'Pu) can be approx- 
imated in trace norm up to an error of 0{h°°) by a semiclassical Fourier integral operator 
with a kernel of the form 



(5.3) 



Here, as in | BK99a |, the phases Si, have to fulfill the Hamilton- Jacobi equations 

K{x,d;,Su{x,^,t)) + dtSu{x,^,t) = 0, Su{x,^,0) = x^. 
The amplitudes a^, G S[!j(l) with asymptotic expansions a^, ~ X^^o determined as 



solutions of certain transport equations 
0{h). Following |pGO0|] further, we choose test functions p E 



i99a\\ with initial conditions a 

'e choc 

supported Fourier transforms p C( 
1 



At=0 ~ X{^u)Pufi + 

with compactly 



oo / 
I 



such that 



h 
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where Tr denotes the operator trace on the Hilbert space 
classical approximation ( ^.3|) one now has to calculate 



Using the semi- 



1 




p{t) J2 tr (5o(x, d,S,)a,,o{x, x, t, ^))ei(^''(-.?.*)-«+^*) dx dt (5.4) 



in leading semiclassical order. This can be done with the method of stationary phase, 
where the stationary points {xi,^st,^u,st,tu,st) of the phase Su{x,^,t) — x^ + Et determine 
periodic points (x^^st, ^i/,st) € ^u,e of the Hamiltonian flow $^ with periods t^^st- Since the 
eigenvalue function X^, is supposed to be non-critical at E, the periods t^^st of the flow $* 
cannot accumulate at zero, see ||Kob87|| . One can hence split p = pi + /)2 in such a way that 
Pi is supported only in a small neighbourhood of zero and p2 = in the vicinity of zero, so 
that the only period in supp pi is the trivial one, ti^^st = 0. The contribution coming from 
Pi to ( p.4| ) is therefore determined by the periodic points with t^^^t = 0. These build up 
the entire level surface Qu,e which, according to assumption (H3,y), is compact. The result 
then reads (see P599| , |EU00l ) 

I 



j 



h 



xiE) 



Pi(0) 

27T 



Oih)). 
(5.5) 



Coming to the contribution of the term with p2 to the expression (|5.4| ), we recall that p2 
has been chosen to vanish in a neighbourhood of zero. The relevant stationary points are 
hence related to periodic orbits of the flow $^ with non-vanishing periods. The condition 
(H6jy) now allows us to employ the methods of [PS99|| , leading to the estimate 



J2xiE,){^P„B^P,)p2 



Ej-E 

h 



(5.6) 



The relations (|5.5| ) and ( |5.6| ) together therefore imply that for every test function p G 
C°°(]R) with Fourier transform p G C^(M) the estimate (|5.5| ) holds with pi replaced by p. 
Hence, the Tauberian argument developed in |[BPU95|| can be applied to yield 



Ej(^I(E,h) 



ijj 



(5.7) 



In this relation one can set the operator B equal to the identity and thus obtains a semi- 
classical expression for the number Ni of eigenvalues of 7i in liE^ K), 



Nr 



i^{E,el{E,h)} = -J2k 

■rr < ^ 



volfi^^s 

-V +o{h 



l-d\ 



'{2TThy 



(5.^ 



where k^, = trPj^^o denotes the dimension of the fibre ranPj^^o = E"^ corresponding to the 
eigenvalue X,y of Hq. The proof is now finished by combining the expressions ( ^.7|) and 

□ 
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Let us add two comments: 



1. Under the additional assumption (Irr^) the Stratonovich-Weyl calculus discussed in 
section || can be applied. It allows to express tr(Pj,o-Bo-Pi/,o) = tr(Vj'i?oV'y) in terms 
of the symbol 60,;^ introduced in (|4.15|) . This then leads to the representation 

volOxJn^^Jo^ (5.9) 

as an integral over the product space Qu,e x Ox. Here the relation k,^ = vol (9a, 
introduced in section ^, enables one to give the right-hand side of (|0| ) a genuinely 
classical interpretation. 



2. The operators B considered in the limit formula ( |5.1| ) have not been restricted to 
those with symbols in the invariant subalgebra S?jy(l) C S° (1). Nevertheless, only 
the diagonal blocks of their principal symbols Bq with respect to the projection 
matrices P,yfi enter on the right-hand side of ( ^.1|) . In particular, this implies that for 
an operator B with a purely off-diagonal principal symbol, i.e., P^flB^P^ Q = for all 
fi = 1, . . . , Z, the semi classical average vanishes. Thus one can replace an operator 
B with symbol B G S°i(l) by its diagonal part '^^VfiBV^, whose symbol is in the 
invariant algebra S?jy(l), without changing the value of the limit on the right-hand 
side of (13). 

So far we have considered expectation values in normalised eigenvectors of Ti. Our intention 
now is to discuss the projections Vyiljj of the eigenvectors of 7i to a fixed almost invariant 
subspace of L'^iW^) ® C". One thus expresses averaged expectation values in the projected 
eigenvectors in terms of classical quantities related to the single Hamiltonian flow In 



order to achieve this one applies Proposition ^71] to operators V^BVy and exploits the 
selfadjointness of Vy. This results in 

Corollary 5.2. Under the assumptions stated in Proposition \5. j| , for each v E {1, . . . , /} 

the restricted limit formula 

teiV; {V.^P„BV.^P,) = . (5.10) 

holds. 

Thus the semiclassical average of the projected eigenvectors Vyipj, with Ej G I{E,h), 
localises on the corresponding level surface flu,E C T*M'^. If one considers ( p.lO| ) for 
different z/, the relative weights of the corresponding projections are determined by the 
relative volumes of the associated level surfaces and the dimensions of the eigenspaces E'^, 
which equal the volumes of the coadjoint orbits Ox- 
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In general, however, the projected eigenvectors Vu'ipj are neither normahsed, nor are 
they genuine eigenvectors of Ti. We therefore now introduce the normahsed vectors 



(5.11) 



Since the projectors Vu only commute with 7i up to a term of 0{h°°), the pairs {Ej, (pj^u) 
are quasimodes with discrepancies rj^^, i.e.. 



This observation only ensures the existence of an eigenvalue of Ti in the interval [Ej — 
rji,,Ej + Tju], which is a trivial statement; it does not imply that (pj^^ is close to an 



eigenvector of 7i, see ||Laz93|| . It therefore is of somewhat more interest to consider the 



operator HVu, whose spectrum inside the interval [E — e,E + e] D I{E, h) is as well purely 
discrete. Following the above reasoning, one then concludes that {Ej,(f)j^t,) is a quasimode 
with discrepancy Vj^^, also for this operator. Thus, if H'Pi/V'jll ^ ch^ for some > and 
hence Vj^^ = 0{h°°), the operator TCVu has an eigenvalue with distance 0{h°°) away from 
Ej. Since there are A^- eigenvalues Ej G I{E, h) one finds as many quasimodes for TiVu- 
But this operator has only 

fc^voU^ 



eigenvalues in I{E,h), compare ( ^.8| ). This observation might suggest that only approx- 
imately Nj of the Ni projected eigenvectors V^i/jj are of considerable size, such that the 
discrepancies of the associated quasimodes are smaller than the distance of Ej to neigh- 
bouring eigenvalues of Ti. This expectation can be strengthened by an application of the 
limit formula ( |5.10| ) with the choice i3 = id. 



which implies that 



Ej&I{E,h) 



One could thus expect that roughly Nj of the projected eigenvectors Vuipj are close to ipj, 
and the rest is such that is semiclassically small. However, ( p.l3| ) does not rule 

out the other extreme situation, provided by projected eigenvectors V^ipj, v = 
equidistributing in the sense that their squared norms are asymptotic to Nj /Nj as — > 0. 
In that case the discrepancies of the associated quasimodes for the operators HVu can be 
estimated as r^^y = 0{h°°). In order now that these quasimodes do not produce more 
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than N'^ eigenvalues of TiVv in /(-E, h), a finite fraction of the eigenvalues Ej of Ti must 
possess spacings to their nearest neighbours of the order ff° . Since in general there exist no 
sufficient lower bounds on eigenvalue spacings, none of the two extreme situations discussed 
above can be excluded so far. 

What is possible, however, is to derive from ( [5 .121 ) an upper bound for the fraction of 
the projected eigenvectors Vyipj that are close in norm to ipj, 

hm G mh)- - ^,11 = o(l)} < ,^-7^^;-^ , 

^-^^ Y.^,=l vol ^^,,E 

see also ||Sch01|| . To obtain lower bounds is notoriously more difficult. The limit formula 
( ^.12| ) only allows to estimate the fraction of projected eigenvectors with norms that tend 
to a finite limit as — > 0. One conveniently measures this fraction in units of the value 
that is expected for equidistributed projections. Therefore, with 5 := ^ ^i^''T'^^7n — ' 
consider 



Since 



Ej(^I{E,h) Ej(^I{E,h) EjeI{E,h) 



|2 



the relative fraction of projected eigenvectors with finite semiclassical limit can be esti- 
mated from below as 

lim^>ilfil^^i^^^. (5.14) 



6 Quantum ergodicity 

Our intention in this section is to consider quantum ergodicity for the normalised eigenvec- 
tors ipj, Ej G I{E, h), of the quantum Hamiltonian Ti. In the case of scalar pseudodiffer- 
ential operators one denotes by quantum ergodicity a weak convergence of the phase space 
lifts of almost all eigenfunctions to Liouville measure on the level surface flE = Hq^{E), 
and proves this to hold if the fiow generated by the principal symbol Hq of the quantum 
Hamiltonian is ergodic on Qe- In the present situation of operators with matrix valued 
symbols, however, each eigenvalue of Hq defines its own classical dynamics. One hence 
can only expect quantum ergodicity to be concerned with statements about the projec- 
tions Vutpj of the eigenvectors to the different almost invariant subspaces of L^(M'^) ® C" 
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(6.1) 



in relation to the behaviour of the associated classical systems. In the preceding section 
we discussed the question of identifying those projected eigenvectors whose norms are not 
semiclassically small. Since presently this problem cannot be resolved directly, quantum 
ergodicity can only be formulated by restricting to those eigenvectors whose squared norms 
exceed a value of S in the semiclassical limit, without specifying them further. 

Conventionally the convergence of quantum states determined by the eigenvectors tpj of 
H is discussed in terms of expectation values of observables in these states. Explicit lifts of 
the eigenfunctions to phase space are then, e.g., provided by their Wigner transforms. The 
choice of the projected eigenvectors V^ipj leads to consider expectation values of diagonal 
blocks VuBVu of operators B with symbols B G S^i(l). On the symbol level the time 
evolution of these blocks is covered by the Egorov theorem p.2| . Representing then the 
blocks of the principal symbols by Stratonovich-Weyl symbols as described in section ^, 
according to Proposition we are faced with the skew-product flows on the product 
phase spaces T*]R'^ x Ox. Since the Stratonovich-Weyl symbols b^^i, defined in equation 
( [4.15| ) that are associated with symbols B G S^j(l) are clearly integrable with respect to 
the measures didi] on the (compact) manifolds Q^^e x 0\, the (assumed) ergodicity of the 
flow implies that 

hm i r{bo,uoY^){x,^,v) dt = —1- / / boAx',e,v') dv'di{x',0 
T^oo T Jo vol Ox Jn,^ J Ox 

holds for almost all initial conditions (x,^,?]) G Qu,e x C^a- In particular, one immediately 
realises that the supposed ergodicity of implies ergodicity for the flow $^ on Q,^ e with 
respect to Liouville measure di. As a consequence the condition (H61,) is automatically 
fulfilled. 

For the subsequent formulation and proof of quantum ergodicity we choose to follow 
in principle the approach of ||Zel96| , |ZZ96|| . This means that we investigate the variance 
of expectation values about their mean in the semiclassical limit. In order to avoid the 
problem of explicitly estimating the norms of projected eigenvectors we here consider the 
normalised vectors (pj^^, defined in (|5.11| ), which have been identified as quasimodes for 
both the operators Ti and TiVu- Moreover, we concentrate on vectors corresponding to 
projected eigenvectors with norms that do not vanish semiclassically, i.e., with > S 

for some fixed 6 G (0, 1). This approach is similar to the one introduced by Schubert 
|Sch01|| in the context of local quantum ergodicity, where an equidistribution was shown for 
quasimodes associated with ergodic components of phase space. In section ^ we estimated 
the relative number N^j/Nj of the associated eigenvectors among all eigenvectors of TC in 
the semiclassical limit from below, see ( ^.14|) . A non-trivial bound could only be obtained 
for 6 < 1 corresponding to 

^ ^ vol Q„^E 
<dy := —I . 

Y.^,=l k,i vol Q^^E 

Therefore, from now on we confine 6 to the interval 6 G (0, 6^), and are thus in a position 
to state our main result. 
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Theorem 6.1. Let Ti be a pseudodifferential operator with hermitian symbol H G S°i(m) 
whose principal part Hq fulfills the conditions (HI) and (H2) of section |^. The eigenvalues 
Ai, . . . , A; of Hq are required to have constant multiplicities and shall obey the conditions 
(H3u)-(H5u) of section^ for all u E {1, . . . , /}. Moreover, they shall be separated according 
to the hyperbolicity condition (HO), 

\X^{x,^) — \fj,{x,^)\ > Cm{x,^) for v ^ ii and + |^| > c. 

Assume now that the symbol H ~ Yl^=o ^"'-^i satisfies the growth condition 

ll^4^)(a^,0llnxn<C„,^ for all e T*R'^ and |«| + |/3|+j> 2 -5,0, dP) 

and that the condition (Irr^) of section holds. If then the flow defined in ^.ISj ) is 
ergodic on Qu,e x C^a with respect to the invariant measure didrj, in every sequence of 
normalised projected eigenvectors {4>j,u}j€N, with WVyipjW^ > 5, 5 E (0,5^,) fixed, one finds 
a subsequence {<t>ja,v}a&i of density one, i.e., 

>^} _^ 

such that for every operator B with symbol B G S°i(l) and principal symbol Bq 

lim(0,„,,, B(j),^,,) = Ms,,,a(&o,.), (6.2) 

where 6o,i/ denotes the Stratonovich-Weyl symbol associated with PufiBoP^yfi. Furthermore, 
the density-one subsequence {(pj^,u}am can be chosen to be independent of the operator B. 

Proof. We start with considering expectation values of the operator B taken in the quasi- 
modes {4>j.u} and denote their variance about the mean ME.u,\{bo,u) of the corresponding 
Stratonovich-Weyl symbol 6o,i/ defined in ( |4.15| ) as 

Sl{E,h):=j^ E |(0.,.,S0.,.)-M^,.,A(V)r- 

'^'^ Ej&I(E,h) 
\\V„i}j\\'^>5 

Due to the definition ( [5 .111) of the normalised vectors (pj^y, this variance can also be written 



as 



Ej£l{E,h) 
||-P.V,IP>'5 
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Allowing for an error of 0{h°°), in this expression the expectation values can be re- 
placed by those of the operator V^iB — MEM,x{bo^u))'Pu whose symbol is in the invari- 
ant subalgebra S['ny(l) C S°i(l). Therefore, since all further requirements are also met, 
the Egorov theorem |3.2| applies and yields that for finite times t G [0, T] the evolution 
U*{t)Vy{B — y^E,v,\{bQ^y))Vi,U{t) of this operator is again a pseudodifferential operator 
with symbol in the class S|?i(l). Taking into account that the ipjS are eigenvectors of Ti 
with eigenvalues Ej, the above expression can be rewritten as 



'^'^ Ej&I{E,h) 



1-2 
3 \\ ' 



where we have defined the auxiliary operator 

B,,T U*{t)Vu{B - ME,uAbo,u))VuKit) dt. (6.3) 

Furthermore, by using using the Cauchy-Schwarz inequality and the lower bound on the 
norms H'Pi/'^/'jll^ > 5 > we obtain as an upper bound 

'^'^ ^ Ej&I(E,h) 

According to equation ( ^.141) the factor Nj/N^j can be estimated from above in the semi- 
classical limit. We hence now consider the semiclassical limit of the expression 

^ Ej£l{E,h) 



to which Proposition 5J. can be applied. To this end one requires the principal symbol 
Bu,T,o of the auxiliary operator B^^t, which follows from Theorem p.2| as 



Bu,T,o = ^ ^ C((P,,o5oP,,o) o dt - M^,,,,(6o,.)P.,o- 

Given this, the limit formula and the estimate (|5.14| ) yield 

1 Yl^^^=lf^^c'^oln^^E vol fi,.,ijtr4,ij (52^0) 



^<5 

I rf , II I *^ — ^ 



lim^2\(^,^) < . ~ 

' {i-5)Kvo\n,,E ELik^yo\n^,E 



^ ^Ms,,,,((symb^^[i?.,T,o])'). 



(6.4) 



when employing the tracial property (v) of Proposition |475 . 
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According to Proposition the Stratonovich-Weyl symbol of B^^t,o can now be easily 
calculated as 

symb^^[S,,T,o(x,0](^) = ^ £ {b^,, oY^){x,^,r]) dt - Me,.^^,^) 



Since we assume the skew-product flow Yj; to be ergodic with respect to didr], the relation 
( |6.1|) implies that symh^^[B,y^T,o{x,^)]{'ri) vanishes in the limit T — oo for almost all points 
{x,^,ri) G ^^^E X Ox- Now, on the right-hand side of (|6.4| ) the square of symb'^^ [i?,^^T,o] 
enters integrated over Qu,e x C>x, so that this expression vanishes as T oo. We hence 
conclude that 



lim SUE, h)=0. 



This, in turn, is equivalent to the existence of a subsequence {4>ja,u}aeN C {0j>}jeN of 
density one, such that equation ( |6.2|) holds. Finally, by a diagonal construction as in 
|Zel87| , |CdV85|] one can extract a subsequence of {(pja,u}aeN that is still of density one in 



{(f)j^i,}j(z-M, such that (|6.2|) holds independently of the operator B. □ 

The version of quantum ergodicity asserted in Theorem |6.1| means that in the semiclassical 
limit the lifts of almost all quasimodes (pj^i, to the phase space T*]R'^ x Ox equidistribute 
in the sense that suitable Wigner functions (weakly) converge to an invariant measure on 
^u,E X C>x that is proportional to didr]. In order to identify the proper Wigner transform 
consider 

= jj ^^^ti[w[c^i^A{^,i)P.{x,i)B{x,i)P.{x,i)) dxdi + 0{K-), 



with the matrix valued Wigner transform 



Wmx,i):-- 



defined for ip G L'^{W'') ® C". We now exploit the Stratonovich-Weyl calculus to conclude 
that on the level of principal symbols 



tr(W^[0i.,.]^.,oSoP.,o) = tr( (KT^oK) 

symb^^[K;W^[0,.,.]K](r7) symb^^[K*i?oK](r7) dr^. 



Ox 



The second factor in the integral has been defined as ho^^ in ( [4.15| ). In analogy to this we 
therefore introduce for ip G L^(]R'^) ® C" the scalar Wigner transform (see also ||BGK01|| ) 



«;,H(x,e,r7) := symb^'^[V;*(x, OW^Ml^^, OK(x, 
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that indeed provides a lift of ip to the phase space T*R x 0\. The statement of Theorem ^ 
can thus be rephrased in that under the given conditions one obtains (in the sense of a 
weak convergence), 

1™ ^^If^Y Wu[<Pu,Ai.^,i,r]) da;d^dr/= ;;;^J^d£(x,0 dr] 

along the subsequence of density one. However, since in (f)j^y the normalisation of Vyipj 
is hidden, an equivalent equidistribution for the lifts of the projected eigenvectors is only 
shown up to a constant. In analogy to the discussion in [^ch01|| this means that in the 



sequence {ipf, Ej G I{E, h)} there exists a subsequence {i/jja} of density one such that as 

{i:,^,V,BV,ij,J = \\V,ij,J\^ME,,,x{bo,u) + o{l), 

with a corresponding statement for the scalar Wigner transforms WylVuipja]- Notice that 
the factor WVui'ja |P is independent of the operator B so that the subsequence can again be 
chosen independently of B. Therefore, a non-vanishing semiclassical limit only exists for 
those subsequences along which the norms H'Py'i/'j^ || do not tend to zero as h 0. These 
subsequences are excluded in the formulation of Theorem |6.1| since 6 is fixed and positive. 

The difficulties with estimating norms of the projected eigenvectors V^ipj arise from 
the presence of several level surfaces ^Iu,e on which the lifts of eigenfunctions potentially 
condense in the semiclassical limit. The situation simplifies considerably, if at the energy 
E all of the / level surfaces except one are empty. 



Corollary 6.2. If under the conditions stated in Theorem \6. j| only the level surface flu,E C 
T*M'^ is non-empty, there exists a subsequence {ipja} of density one in {ipj; Ej G I{E,h)}, 
independent of the operator B, such that 

\im{ipj^,V^BVi,tpjJ = S^uME.uAbo^u)- 

In this situation the norms HP^^Aj^:, || converge to one for /i = z/ and to zero otherwise as 
h —>■ along the subsequence. The lifts of the eigenvectors therefore condense on the only 
available level surface in T*R'^, as one clearly would have expected. 

Remark 6.3. As a condition for quantum ergodicity to hold we have assumed the skew- 
product fiow Yj; on Q^^e x C'a to be ergodic. The reason for introducing this fiow was to 
formulate a genuinely classical criterion in terms of a dynamics on the symplectic phase 
space T*]R'^ x Ox. The formulation will be somewhat simpler, if one refrains from insisting 
on a completely classical description and employs the skew-product fiow Yj; defined on 
T*]R'^ X G, see ( |4.6| ), instead. Then the use of the Stratonovich-Weyl calculus can be 
avoided. Such a formulation is based on a hybrid of the classical Hamiltonian fiow on 
T*]R^ and the dynamics represented by the conjugation with the unitary matrices D^,, which 
appears to be quantum mechanical in nature. Both formulations, however, are equivalent 
in the sense that, first, the Stratonovich-Weyl calculus relates the quantum dynamics in the 
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eigenspace to a classical dynamics on the coadjoint orbit in a one-to-one manner. Second, 
in appendix^ we show that the skew-product Y^; on VL^e x 0\ is ergodic, if and only if the 
skew-product Y^, is ergodic on VL^^e x G. One can therefore formulate Theorem without 
recourse to the Stratonovich-Weyl calculus once the limit y^E,u,\{bQ^y) is expressed as 

y^E,u^{ho,u) = ^tr4,£(P.,o5oP.,o), 

K,, 



see ( ^.91) . Up to equation (|6.4|) the proof of Theorem |6.1| proceeds in the same manner as 
shown. From this point on one can then basically follow the method of [ [BGOO| | , and to this 



end represents the principal symbol -Bi/,t,o of the auxiliary operator ( |6.3|) in terms of the 



isometrics V^, 



V:B,,t,oV, = ^ r D:{{V:BoV,) o $t)Z}, dt - ltr4,i=;(K:5oK). 

-L Jo f^u 

We now suppose that the flow F^f is ergodic on Qi^^e x G and choose the function F{x, ^,g) := 
p{gy{V;BoV^){x,^)p{g) E L^{Vt^^E x G) ®Mfc,^(C) to exploit the ergodicity. This yields for 
almost all initial values (x, ^, g) G VL^^ e x G that 

hm p{grv:{x, OBu,tM Op{9) 

I — >oo 

p{hnv:BoV,){y,()p{h) dhdiivX) -^tTi,,E{V:BoK). 

Furthermore, since the representation (p, C'^") is assumed to be irreducible and the integral 
in the above expression is invariant under conjugation with arbitrary elements of U{k^), 
Schur's lemma implies that this integral is a multiple of the identity in C'^", leading to 

pihnv:BoV,){y,C)p{h) dhde{y,() = ^ tr 4,i^(K*5oK). 



Due to the way the principal symbol -Bi/,t,o enters on the right-hand side of ( |6.4]) , the 
conjugation with Vu{x,^)p{g) as well as the restriction to almost all {x,^,g) is inessential, 
so that again one concludes a vanishing of S2^{E, h) as h ^ 0. 
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Appendices 



A Relations for Poisson brackets of matrix valued 
functions 

In this appendix we collect some relations for Poisson brackets of matrix valued functions 
on the phase space T*R'^ that are needed in section |[ These relations are already stated 
in |[EW96| , pMMP9"7| , |SpoOCI|| and can be verified by straightforward calculations. 

Our convention for the Poisson bracket of smooth matrix valued functions A, B ^ 

C°°(T*M'^) (g) M„(C) is 

{A, B} := d^A d^B - d^A d^B. 

The first general relation then reads 

A{B, C} - {A, B}C = {AB, C} - {A, BC}. (A.l) 

Furthermore, for the projection matrices P = PP one finds 

P{\,P}P = 0, (A.2) 

where A is any smooth scalar function on T*M'^. 

For B G C°°(T*M'^) (g) M„(C) commuting with P one then derives 

P{X, B}P = {A, PBP} - [PBP, [P, {A, P}]], (A.3) 

In particular, using ([A. 1|) for projection matrices one obtains 

P{P, B} - {P, P}B = {P, B} - {P, PB} 

and 

B{P, P} - {B, P}P = {BP, P} - {B, P}. 

Using these relations together with the condition [B,P] = one gets 

P{{B, P} - {P, B})P = [B, P{P, P}P]. 

Furthermore, for different projection matrices and P^, with P^Pi^ = for u fi the 
general relation ( [A. 1| ) implies 

P^,{Pu,Pu} = -{P„Pu}{l-Pu) 

and 

{P„PJP^ = -(1-P,){P„P^}. 
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In the case [P^, -B] = = [P^, B] one finds 

P,{P,„ B} - {P,, P^}B = -{P., P^fi}, 
P{P^, P,} - {P, P^}P, = {PP^, P,}. 

These equations imply 

P,({P,P^} - {P,„P})P. = -[P,P4P^,P4PJ. 
One can now apply the above relations to expressions of the type arising in section ^, i.e. 

P^ [j^Bit) + i({p(t), A.PJ - {A,P,,P(t)}) +i[P(t),Pi]^ P^ 
= §-^P,B{t)P, - 6,,{X,, P,B{t)P^} 



+ 



y (-1)^-P4P,, P,}P,, - S,^[P,, {A,, PJ] - iP^i/lP^, P^P(t)P^ 



Therefore, the definition 



Pi := i(-l)'^-^P^{P,, P,}P, - iS,,[P,, {A,, P,}] + P^PiP^ (A.4) 
allows to conclude that 

^P^PP^ - 6,^{X,, P^BP^} - i[Pi, P^PP^] = 0. (A.5) 

B A relation between the ergodicity of two skew-pro- 
duct flows 

In section ^ we considered two types of skew-product dynamics built over the Hamiltonian 
flows $^ on T*M'^. Both derive from the dynamics in the eigenvector bundles E'^ — > T*R'^ 
given by conjugating the diagonal blocks of principal symbols with the transport matrices 
di^u along integral curves of the Hamiltonian flows. After having fixed local orthonormal 
bases in the fibres, or isometries V,^{x,^) : C'^" E'^{x,^), respectively, the transport 
matrices d^^ have been represented by unitary ki, x ki, matrices D^, leading to the skew- 
product flows on T*M'^ x \J{ki,). We then noticed that the dynamics in the fibres might 
not exhaust the whole group \J{k^), but only some subgroup G, which is then represented 
in U(fc^). This led us to consider the skew-product flows on T*M^ X G, given as 
Y^{x,C,,g) = {^l,{x,^), gu{x,^,t)g), see (|4.5| ) and ( [4.61) . Assuming that the representation 



p of G in U(/ci/) is irreducible, we constructed a representation of the flbre dynamics on 
the coadjoint orbit Ox of G determined by p. We thus arrived at the skew-product flows 
Fj on the symplectic phase spaces T*R'^ x Ox, with Y^{x,^,ri) = (<l>^(a;, ,^), Ad*^(^.^g r/), 
see ( [4.12[ ) and ( |4.13| ). In section ^ we required either the flows Y^ or yj, restricted to the 
level surfaces ^Iu,e C T*R'^ in the base manifold, to be ergodic relative to the respective 
invariant measures didg or d£d?7. We now show: 
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Proposition B.l. The flow "Kj : Q,„^e x G - 
and only if the associated flow Y^; : Qu,e x C^p 



Qu,E X G is ergodic with respect to didg, if 
Qu,E X Oi, is ergodic with respect to didf]. 



Proof. A convenient characterisation for tlie ergodicity of a flow $* on a probability space 
(S, dm) with invariant measure dm employs the flow-invariant subsets of S: The flow is 
ergodic with respect to dm, if and only if every measurable flow-invariant set has either 
measure zero or full measure. We now first consider the 'if direction asserted in the 
proposition and to this end assume that on Q^^; x 0\ is ergodic with respect to didr]. 
Hence every measurable Fj-invariant set B C ^Iu,e x Ox has either measure zero or full 
measure. In order to relate these sets with subsets of Qu,e x G we recall the composed 
map G G/G\ 0\ from section ^, where vr denotes the canonical projection of G onto 
G/Gx and k is the diffeomorphism that identifies G/G\ with 0\. One then realises that 
the following diagram commutes: 



{Kix,0,9u{x,tt)g) 



9uix,^,t)gGx) 



id„ 



iB.l] 



where is induced by under idT-md xtt. According to this diagram a yj-invariant 
set A C Qu,E X G projects to a F^^-invariant subset (idT.Rd xk o n){A) of Q^^e x Ox. The 
assumed ergodicity of then implies that the measure of (idx*Md x/s;o7r)(y4) is zero or one. 
Now the normalised Haar measure dg on G projects under k o tt to the volume measure drj 
on the coadjoint orbit Ox- This can be obtained from the Fubini theorem (cf. | [BtD85| |) 
which states for every / G L^{Ox) that 



7T*K*f)ig) dg= / iK*f) o nigh) dh ] d{gGx) 

G JG/Gx \JGx 



(B.2) 



{n*f){gGx) d{gGx). 



G/Gx 



Here dh denotes the normalised Haar measure on Gx and d{gGx) is the normalised left 
invariant volume form on G/Gx arising from the volume form on the coadjoint orbit under 
the pullback k*. Hence, the sets A and (idT*Rti XK,o'n-)[A) have identical measures and thus 
the measure of A is either zero or one. Therefore, the assumed ergodicity of Yj; implies 
ergodicity of Y^. 

In order to prove the opposite direction one simply reverses the above argument: Start- 
ing with yj-invariant subsets of Qu,e x Ox, one lifts these to Qu,e x G. Due to the commuting 
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diagram ( [B.l| ) these lifts are yj-invariant and therefore, according to the assumed ergodic- 
ity of yj, have measure zero or one. Again the Fubini theorem (p.2| ) imphes equal measures 
of the sets and their lifts. Hence is ergodic. □ 
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